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Abstract 

The problem of optimal resource allocation is studied for ergodic fading orthogonal multiaccess 
relay channels (MARCs) in which the users (sources) communicate with a destination with the aid of a 
half-duplex relay that transmits on a channel orthogonal to that used by the transmitting sources. Under 
the assumption that the instantaneous fading state information is available at all nodes, the maximum 
sum-rate and the optimal user and relay power allocations (policies) are developed for a decode-and- 
forward (DF) relay. With the observation that a DF relay results in two multiaccess channels, one at 
the relay and the other at the destination, a single known lemma on the sum-rate of two intersecting 
polymatroids is used to determine the DF sum-rate and the optimal user and relay policies. The lemma 
also enables a broad topological classification of fading MARCs into one of three types. The first type 
is the set of partially clustered MARCs where a user is clustered either with the relay or with the 
destination such that the users waterfill on their bottle-neck links to the distant receiver The second 
type is the set of clustered MARCs where all users are either proximal to the relay or to the destination 
such that opportunistic multiuser scheduling to one of the receivers is optimal. The third type consists 
of arbitrarily clustered MARCs which are a combination of the first two types, and for this type it is 
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shown that the optimal poHcies are opportunistic non-waterfiUing solutions. The analysis is extended to 
develop the rate region of a i^T-user orthogonal half-duplex MARC. Finally, cutset outer bounds are used 
to show that DF achieves the capacity region for a class of clustered orthogonal half-duplex MARCs. 

Index Terms 

Multiple-access relay channel (MARC), decode-and-forward, ergodic capacity. 

I. Introduction 

Node cooperation in multi-terminal wireless networks has been shown to improve performance by 
providing increased robustness to channel variations and by enabling energy savings (see [l]-[7] and the 
references therein). A specific example of relay cooperation in multi-terminal networks is the multi-access 
relay channel (MARC). The MARC is a network in which several users (source nodes) communicate with 
a single destination with the aid of a relay [8]. The coding strategies developed for the relay channel [9] 
extend readily to the MARC [10]. For example, the strategy of [9, Theorem 1], now often called decode- 
and-forward (DF), has a relay that decodes user messages before forwarding them to the destination [3], 
[11]. Similarly, the strategy in [9, Theorem 6], now often called compress-and-forward (CF), has the 
relay quantize its output symbols and transmit the resulting quantized bits to the destination [10]. 

We consider a MARC with a half-duplex wireless relay that transmits and receives in orthogonal 
channels. Specifically, we model a MARC with a half-duplex relay as an orthogonal MARC in which 
the relay receives on a channel over which all the sources transmit, and transmits to the destination on 
an orthogonal channel^. This channel models a relay-inclusive uplink in a variety of networks such as 
wireless LAN, cellular, and sensor networks. The study of wireless relay channels and networks has 
focused on several performance aspects, including capacity [1], [3], [9], diversity [2], [4], [13], outage 
[14]-[16], and cooperative coding [17], [18]. Equally pertinent is the problem of resource allocation in 
fading wireless channels where both source and relay nodes can allocate their transmit power to enhance 
a desired performance metric when the fading state information is available. Resource allocation for a 

'Yet another class of orthogonal single-source half-duplex relay channels is defined in [12] where the source and relay transmit 
on orthogonal bands. The source transmits in both bands, one of which is received at the relay and the other is received at the 
destination, such that the relay also transmits on the band received at the destination. In contrast to [12], we assume that all 
sources transmit in only one of orthogonal bands and the relay transmits in the other. Furthermore, we assume that signals in 
both bands are received at the destination. Later in the sequel we briefly discuss the general model where the sources transmit 
on both bands. 
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variety of relay channels and networks has been studied in several papers, including [5], [14], [19]-[21]. 
A common assumption in all these papers is that the source and relay nodes are subject to a total power 
constraint. 

For a wireless fading relay chaimel, i.e., a single-user specialization of a fading MARC, the problem 
of resource allocation when the source and relay nodes are subject to individual power constraints in 
studied in [6] (see also [22]). The authors formulate the problem as a max-min optimization. They draw 
parallels with the classical miiumax optimization in hypothesis testing to show that, depending on the 
joint fading statistics, the resource allocation problem results in one of three solutions. The three solutions 
broadly correspond to three types of channel topologies, namely, source-relay clustering, relay-destination 
clustering, and the non-clustered (arbitrary) topology. 

Resource allocation in multiuser relay networks has been studied recently in [23]-[25]. The authors in 
[23] and [25] consider a specific orthogonal model where the sources time-duplex their transmissions and 
are aided in their transmissions by a half-duplex relay, while in [24] the optimal multiuser scheduling is 
determined under the assumption of a non-fading backhaul chaimel between the relay and destination. 
In contrast, in this paper, we consider a more general multiaccess chaimel with a half-duplex relay and 
model all inter-node wireless links as ergodic fading channels with perfect fading information available at 
all nodes. Assuming a DF relay, we develop the optimal source and relay power allocations and present 
the conditions under which opportunistic time-duplexing of the users is optimal. 

The orthogonal MARC is a multiaccess generalization of the orthogonal relay channel studied in [6] ; 
however, the optimal DF policies developed in [6] do not extend readily to maximize the DF sum-rate 
of the MARC. This is because unhke the single-user case, in order to determine the DF sum-rate for the 
MARC, we need to consider the intersection of the two multiaccess rate regions that result from decoding 
at both the relay and the destination. Here, we exploit the polymatroid properties of these multiaccess 
regions and use a single known lemma on the sum-rate of two intersecting polymatroids [26, chap. 46] 
to develop inner (DF) and outer bounds on the sum-rate and the rate region and specify the sub-class of 
orthogonal MARCs for which the DF bounds are tight. 

A lemma in [26, chap. 46] enables us to classify polymatroid intersections broadly into two sets, 
namely, the sets of active and inactive cases. An active or an inactive case result when, in the region of 
intersection, the constraints on the i^T-user sum-rate at both receivers are active or inactive, respectively. 
In the sequel we show that inactive cases suggest partially clustered topologies where a subset of users is 
clustered closer to one of the receivers while the complementary subset is closer to the remaining receiver. 
On the other hand, active cases can result from specific clustered topologies such as those in which all 
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sources and the relay are clustered or those in which the relay and the destination are clustered, or more 
generally, from arbitrarily clustered topologies that are either a combination of the two clustered models 
or of a clustered and a partially clustered model. For both the active and inactive cases, the polymatroid 
intersection lemma yields closed form expressions for the sum-rates which in turn allow one to develop 
the sum-rate optimal power allocations (policies). 

We first study the two-user orthogonal MARC and develop the DF sum-rate maximizing power policies. 
Using the polymatroid intersection lemma we show that the fading-averaged DF sum-rate is achieved 
by either one of five disjoint cases, two inactive and three active, or by a boundary case that lies at the 
boundary of an active and an inactive case. We develop the sum-rate for all cases and show that the sum- 
rate maximizing DF power policy either: 1) exploits the multiuser fading diversity to opportunistically 
schedule users analogously to the fading MAC [27], [28] though the optimal multiuser policies are not 
necessarily water-filling solutions, or 2) involves simultaneous water-filling over two independent point- 
to-point links. Using similar techniques, we also develop the two-user DF rate region. 

Next, we generalize the two-user sum-rate and rate region analysis to the K-usei channel and show 
that the inactive, active, and boundary cases correspond to partially clustered, clustered, and arbitrarily 
clustered topologies, respectively. Finally, we develop the cutset outer bounds on the sum-capacity of an 
ergodic fading orthogonal and non-orthogonal i^-user Gaussian MARC. We show that DF achieves the 
sum-capacity for a class of half-duplex MARCs in which the sources and relay are clustered such that 
the outer bound on the i^-user sum-rate at the destination dominates all other sum-rate outer bounds. 
We also show that DF achieves the capacity region when the cutset bounds at the destination are the 
dominant bounds for all rate points on the boundary of the outer bound rate region. 

In the course of developing the main results of this paper, we also show that DF achieves the capacity 
region of a class of degraded discrete memoryless and Gaussian non-fading orthogonal MARCs where 
the received signal at the destination is physically degraded with respect to that at the relay conditioned 
on the transmit signal at the relay. The relatively few capacity results known for specific classes of 
full-duplex single-user relay channels, such as those for degraded relay channels [9, Theorem 5] and for 
a class of orthogonal relay channels [12], have not been straightforward to extend to the MARC. The 
result developed here is the first in which the entire capacity region is given for a class of degraded 
MARCs. In contrast, in [29] it is shown that DF achieves the sum-capacity of a class of full-duplex 
degraded Gaussian MARCs for which the polymatroid intersections at the relay and destination belong 
to the active set. 

The paper is organized as follows. In Section |lll we present the channel models and introduce 
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polymatroids and a lemma on their intersections. In Section |nl] we develop the DF rate region for 
ergodic fading orthogonal MARCs. In Section JV] we develop the power policies that maximize the DF 
sum-rate for a two-user MARC. We extend the analysis to the K-user orthogonal MARC as well as to 
non-orthogonal models in Section |Vl In Section IVIII we present outer bounds and illustrate our results 
numerically. We summarize our contributions in Section IVIIII 

II. Channel Model and Preliminaries 

A. Orthogonal Half-Duplex MARC 

A K-usei MARC consists of K source nodes numbered 1,2, . . . , K, a relay node r , and a destination 
node d. We write /C = {1,2,..., K} to denote the set of sources, T = /C U {r} to denote the set of 
transmitters, and V = {r, d} to denote the set of receivers. In an orthogonal MARC, the sources transmit 
to the relay and destination on one channel, say channel 1, while the half-duplex relay transmits to the 
destination on an orthogonal channel 2 as shown in Fig. [T] Thus, a fraction 9 of the total bandwidth 
resource is allocated to channel 1 while the remaining fraction 9 = 1 — 9 is allocated to channel 2. In 
the fraction 9, the source k, for all k G )C, transmits the signal while the relay and the destination 
receive Yr and i respectively. In the fraction 9, the relay transmits Xr and the destination receives 
Yii,2 where the sources precede the relay in the transmission order. In each symbol time (channel use), 
we thus have 

= Ylk=l Hr,kXk + Zr, (1) 

Yd,i = ELiHd,kXk + Zd,i, and (2) 

where Zr, Z^ i, and Zd 2 are independent circularly symmetric complex Gaussian noise random variables 
with zero means and unit variances. We write H_ to denote a vector of fading gains, Hk^m, for all k G 
V and m G T, k ^ m, such that his a realization for a given channel use of a jointly stationary and 
ergodic (not necessarily Gaussian) fading process {H}. Note that the channel gains H^^m^ for all k, m, 
are not assumed to be independent. We assume that the fraction 9 is fixed a priori and is known at all 
nodes. Since the relay is assumed to be causal, we note that the signal Xr at the relay in each channel 
use depends causally only on the Yr received in the previous channel uses. 

Over n uses of the channel, the source and relay transmit sequences {X^^i} and {Xr^i}, respectively. 
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are constrained in power according to 

n 

\Xk/ < nPk , for all keT. (4) 

1=1 

Since the sources and relay know the fading states of the links on which they transmit, they can allocate 
their transmitted signal power according to the channel state information. We write Pk{H_) to denote 
the power allocated at the fc*^ transmitter, for all A; € T, as a function of the channel states H_. For an 
ergodic fading channel, ^ then simplifies to 

IE [Pk{H)] < Pk for all k£T (5) 

where the expectation in (O is over the distribution of H_. We write P {H_) to denote a vector of power 
allocations with entries P^ {H_) for all k € T, and define V to be the set of all P {H_) whose entries 
satisfy (HI). Throughout the sequel, we refer to the fractions 6 and 1 — ^ as the first and second fractions, 
respectively. 

B. Polymatroids 

In the sequel, we use the properties of polymatroids to develop the ergodic sum-rate results. Poly- 
matroids have been used to develop capacity characterizations for a variety of multiple-access channel 
models including the MARC (see for e.g., [11], [28], [30]). Furthermore, in [30], Han demonstrates 
that for certain multi-terminal channels, polymatroid intersections need to be considered. To the best 
of our knowledge, this is the first work where the polymatroid intersection lemma has been used to 
explicitly characterize sum-rates and sum-capacity, where possible. We review the following definition 
of a polymatroid. 
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Definition 1: Let /C = {1, 2, . . . , K} and / = 2^ — > be a set function. The polyhedron 

B{f) = {{Ri,R2,...,RK):Rs<f{S), forall5C/C, Rk>0} (6) 

is a polymatroid if / satisfies 

1) / (0) = (normalization) 

2) / (5) < / (P) if 5 C P (monotonicity) 

3) / (5) + / (P) > / (5 U 7?) + / (5 n V) (submodularity). 

Remark 1: The submodularity property in Definition [T] above is equivalent to requiring, for all ki,k2 
in /C with ki ^ k2, h ^ S, k2 ^ S, that / satisfies [26, Ch. 44] 

fiS U {A:i}) + fiS U {^2}) > f{S) + fiS U {ki, k2}. (7) 

This property is used in [11] to show that the rate regions achieved at both the relay and the destination 
in a full-duplex MARC are polymatroids. 

We use the following lemma on polymatroid intersections to develop optimal inner and outer bounds 
on the sum-rate for K-user half-duplex MARCs. 

Lemma 1 ( [26, p. 796, Cor. 46.1c]): Let Rs < fi (S) and Rs < f2 (S), for all 5 C /C, be two 
polymatroids such that /i and /2 are nondecreasing submodular set functions on /C with fi (0) = /2 (0) = 
0. Then 

max Rk = min (/i (5) + ^ (/C\5) ) . (8) 
SQJC 

Lemma \T\ states that the maximum i^-user sum-rate Rjc that results from the intersection of two 
polymatroids, Rs < /i (5) and Rg < f2 (5), is given by the minimum of the two ET-user sum-rate 
planes /i (/C) and /2 (/C) only if both sum-rates are at most as large as the sum of the orthogonal rate 
planes /i (5) and /2 (/C\5), for all 7^ 5 C /C. We refer to the resulting intersection as belonging to 
the set of active cases. 

When there exists at least one 7^ 5 C /C for which the above condition is not true, an inactive case 
is said to result. For such cases, the maximum sum-rate in ([S]) is the sum of two orthogonal rate planes 
achieved by two complementary subsets of users. As a result, the i^-user sum-rate bounds /i(/C) and 
/2(/C) are no longer active for this case, and thus, the region of intersection is no longer a polymatroid 
with 2^ — 1 faces. For a K-user MARC, there are 2^ — 2 possible inactive cases. 

The intersection of two polymatroids can also result in a boundary case when for any 5 C /C, /i (S) + 
/2 {IC\S) is equal to one or both of the K-user sum-rate planes. The orthogonality of the planes /i (5) 
and /2 {IC\S) implies that no two inactive cases have a boundary and thus a boundary case always arises 
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between an inactive and an active case. Note that by definition, a boundary case is also an active case 
though for ease of exposition, throughout the sequel we explicitly distinguish between them. From ([D, 
there are three possible active cases corresponding to the three cases in which the sum-rate plane at one 
of the receivers is smaller than, larger than, or equal to that at the other. In fact, the case in which the 
sum-rates are equal is also a boundary case between the other two active cases. Thus, there are a total 
of (2^ — 1) boundary cases for each active case. 

In summary, the inactive set consists of all intersections for which the constraints on the two sum-rates 
are not active, i.e., no rate tuple on the sum-rate plane achieved at one of the receivers lies within or on 
the boundary of the rate region achieved at the other receiver. On the other hand, the intersections for 
which there exists at least one such rate tuple such that the two sum-rates constraints are active belong to 
the active set. Thus, by definition, the active set also includes those boundary cases between the active 
and inactive cases for which there is exactly one such rate pair. 

III. Two-User Orthogonal MARC: Ergodic DF Rate Region 

The DF rate region for a discrete memory less MARC and a full-duplex relay is developed in [3, 
Appendix A] (see [11] for a detailed proof). For this model, X^, k £ T, denotes the transmit signals at the 
sources and relay and Yr and Yd, denote the received signals at the relay and destination, respectively. The 
rate region is achieved using block Markov encoding and backward decoding. The following proposition 
summarizes the DF rate region. 

Proposition! ( [11, Appendix I]): The DF rate region is the union of the set of rate tuples i?2) ••• j 
Rk) that satisfy, for all 5 C /C, 

Rs < min {I{Xs; Yr\Xs^VicXrU), I{XsXr; Yd\Xs^Vs^U)} (9) 

where the union is over all distributions that factor as 

vi"^) ■ (llk=iPMu)p{xk\vk,u)^ -pixrlvicu) ■p{yr,yd\xr)- (10) 
Remark 2: The time-sharing random variable U ensures that the region of Theorem [J is convex. 
Remark 3: The independent auxiliary random variables Vk, k = 1,2, . . . , K, help the sources cooperate 
with the relay. 

In [10] (see also [31, Proposition 2.5]), the DF rate bounds for a discrete memoryless MARC with a 
half-duplex relay are developed. For the orthogonal MARC model studied, since the sources and relay 
transmit on orthogonal channels, the need for auxiliary random variables Vk, for all k, that model the 
coherent combining gains is eliminated. Under the assumption that the transmit (bandwidth) fractions 6 
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and 1 — 6 at the users and relay, respectively, are known at all nodes, the following proposition summarizes 
the DF rate region for the orthogonal half-duplex MARC. 

Proposition 2: The DF rate region of a orthogonal MARC is the union of the set of rate tuples 
i?2, • • • , Rk) that satisfy, for all 5 C /C, 

Rs < mm{ei{Xs;Yr\Xs^U),ei{Xs;Yd,i\Xs^U) + 9I{Xr;Yd,2\U)} (11) 

where the union is over all distributions that factor as 

piu) ■ (^9 Y\^^-^p{xk\u) ■ piyr,yd\xic) + Op{xr\u) ■ p{yd\xr)^ . (12) 
Definition 2: A parallel MARC is a collection of M MARCs, for which the inputs and outputs of 

parallel channel (sub-channel) j, j = 1,2, . . . ,M, are X^j, € /C U {r} and (Yrj,Ydj), respectively, 

such that conditioned on its inputs, the outputs of each sub-channel are independent of the inputs and 

outputs of other sub-channels. 

Theorem 1: For the parallel MARC, the DF rate region is the union of the set of rate tuples i?2, • • • , 

Rk) that satisfy, for all 5 C /C, 

Rs < min < I{Xs,m', Yr^ml^S" ,mV}C,mXr,mUm) , 

EM ^ 
I{-^S.,mXr,m',Ydm\XS'=,myS'',mUm) } (13) 
m=l J 

where the union is over all distributions that factor as 

) ■ (14) 

Proof: The inner bounds in ([T3]) are obtained by setting U = {Ui U2 ... UM),Vk = (Vk^iVk^2 ••• V/c^m), 
Xk = {Xk,i Xk,2 ■ ■ ■ Xk^u), Yr = {Yr^i Yr^2 ■ ■ ■ 5^r,M), and Yd = {YdA Ydfl . . . YdM), in © and 
choosing {Um, VjCm, Xjc^m) to be independent for all m. ■ 
For the (half-duplex) orthogonal Gaussian MARC with a fixed H_ and 9 that is assumed known at 
all nodes, we consider Gaussian signaling with zero mean and variance at transmitter k such that 
Xk ~ CM{0,Pk), for all k ^ K,. Thus, from (fTTl l the DF rate region includes the set of all rate pairs 
(fii,i?2) that satisfy 

R, < min (^I^^lpi^ + oc (B^L^^ , OC ^Il^^A^ | , , = 1, 2 (15) 

and 

R^ + R2< min !^9C (j^ ^^^^^ + 9C (^I^h^^ , ec (^g L^-^^ | . (I6) 
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For a stationary and ergodic process {H}, the channel in ([T])-© can be modeled as a set of parallel 
Gaussian orthogonal MARCs, one for each fading instantiation H_. For a power policy P_{H), the DF 
rate bounds for this ergodic fading channel are obtained from Theorem [J by averaging the bounds in 
([13] ) and ([T6l ) over all channel realizations. The ergodic fading DF rate region, TZdf, achieved over all 
P_{H_) € V, for a fixed bandwidth fraction 9, is summarized by the following theorem. 

Theorem 2: The DF rate region TZdf of an ergodic fading orthogonal Gaussian MARC is 



7^DF= U {T^r{P)r\nd{P)} 



(17) 



where, for all 5 C /C, we have 



UriP) 



and 



(i?i,i?2) ■■Rs<K 



9C 



( E \Hr,k\ 
k&S 



PkiK) 



\ 



(18) 



T^d (P) 



{Ri,R2) : Rs <^ 



ec 



/E \Hd,k?Pkm\ 

k£S 



\ 



+ 9C 



\HdA Pr{H) 



J 



(19) 



Proof: The proof follows from the observation that the channel in ([TJ-© can be modeled as a set 
of parallel Gaussian orthogonal MARCs, one for each fading instantiation H_. Thus, from Theorem [T] for 
Gaussian inputs and for each P{H_) € V, the regions TZr {P.{H.) ) and TZd {P{H_) ) are given by the bounds 
in ([TSl l and ( fT9l ). respectively. The DF rate region, IZdf, is given by the union of such intersections, one 
for each P,{H_) G V. The convexity of TZdf follows from the convexity of the set V and the concavity 
of the log function. Consider two rate tuples {R'l, R'2) and (i?'/, R2) that result from the policies P_'{H_) 
and P_"{H_), respectively. For any A > such that A = 1 — A, and for all A; = 1, 2, from (fTSl) . we bound 



Rk = XR'j^ + (l — X) R'l achieved at the relay as 



< OE 



Hr,k\^ Pk{H) \ 



< OE 




+ xeE 



c 



xpm)+^Pkm) 



(20) 
(21) 
(22) 



where ((2TI) follows from Jensen's inequality and (|22]) follows from the convexity of the set V such that 
P(H) = (XP'(H) + XP"{H_)) € v. Thus, we see that the bound on R^ is achievable. One can similarly 
bound the sum-rate Ri + R2 achieved at the relay thus proving that the tuple {Ri,R2) € TZ-r. The same 
approach also allows us to show that R2) € TZd, thus proving that Rdf is convex. ■ 
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Proposition 3: TZr (P(H)) and TZd (P(H)) are polymatroids. 

Proof: In [11, Sec. IV.B], it is shown that for each choice of the input distribution in (fTOl) . the DF 
rate region in ^ is an intersection of two polymatroids, one resulting from the bounds at the relay and 
the other from the bounds at the destination. For the orthogonal MARC, the bounds in ([TT]) . relative to 
(|9l), involve a weighted sum of mutual information expressions; using the same approach as in [11, Sec. 
IV.B], the submodularity of these expressions can be verified in a straightforward manner. ■ 

Remark 4: The DF rate region given by (fTSl) and ([T6l) is achieved using block Markov encoding at 
the sources. For the ergodic fading model, the rates in Theorem |2] are obtained assuming that all fading 
instantiations are seen in each such block. 

In the following section, we develop sum-rate optimal DF power policies. 

IV. Two-User Orthogonal MARC: DF Sum-Rate Optimal Power Policy 

For ease of notation, throughout the sequel, we write Rj(j to denote the sum-rate bound on the users in 
A and 7^^™ to denote the sum-rate obtained by successively decoding the users in IC\A before decoding 
those in A at receiver j = r, d. For the two-user case, Rjcj and Raj^ for all ^ C /C are given by the 
sum-rate and single-user bounds in ( [T8] ) and ( [T9l ) at the relay and destination, respectively. The rate TZj^^'j 
= Ric,j — Ric\A,j> for all ^ C /C, is obtained by successively decoding the users in ]C\A before decoding 
those in A at the corner points of the regions TZr and TZa- 

The region TZdf in ( flTl ) is a union of the intersections of the regions 'R,r{P_{lP)) and TZd(P(H)) 
achieved at the relay and destination respectively, where the union is over all P_{H) € V. Since TZdf 
is convex, each point on the boundary of TZdf is obtained by maximizing the weighted sum jiiRi + 
/i2^2 over all P.{H_) € 'P, and for all /xi > 0, > 0. Specifically, we determine the optimal policy 
P_*{H_) that maximizes the sum-rate Ri + R2 when fii = ^2 = 1- Observe from (ITtI) that every point 
on the boundary of TZdf results from the intersection of the polymatroids (pentagons) TZr{P_{H_)) and 
TZd{P{H_)) for some P(H). In Figs. |2] and [3] we illustrate the five possible choices for the sum-rate 
resulting from such an intersection for a two-user MARC of which two belong to the inactive set and 
three to the active set. 

The inactive set consists of cases 1 and 2 in which user 1 achieves a significantly larger rate at the 
relay and destination, respectively, than it does at the other receiver; and vice-versa for user 2. The active 
set includes cases 3a, 36, and 3c shown in Fig. |2] in which the sum-rate at relay r is smaller, larger, or 
equal, respectively, to that achieved at the destination d. The three boundary cases between case 1 and 
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Fig. 2. Rate regions Rr{P) and Rd{P) and sum-rate for case 1 and case 2. 



the three active cases are shown in Fig. |4] while the remaining three between case 2 and the active cases 
are shown in Fig. |5] We denote a boundary case as case (/, n), / = 1, 2, n = 3a, 36, 3c. 

We write Bi V and Bi^n Q V to denote the set of power poUcies that achieve case i, i = 
1, 2, 3a, 36, 3c, and case {l,n), I = 1,2, n = 3a, 36, 3c, respectively. We show in the sequel that the 
optimization is simplified when the conditions for each case are defined such that the sets Bi and Bi^n 
are disjoint for all i, I, and n, and thus, are either open or half-open sets such that no two sets share a 
boundary. Observe that cases 1 and 2 do not share a boundary since such a transition (see Fig. O requires 
passing through case 3a or 36 or 3c. Finally, note that Fig. [3] illustrates two specific TZr and TZd regions 
for 3a, 36, and 3c. For ease of exposition, we write B^ = B^a U B-^h U B^c, where Bi, i = 3a, 36, 3c. 

In general, the occurrence of any one of the disjoint cases depends on both the channel statistics and 
the policy P,{H_). Since it is not straightforward to know a priori the power allocations that achieve a 
certain case, we maximize the sum-capacity for each case over all allocations in P and write P_^'^\H_) 
and P^''''^\H_) to denote the optimal solution for case i and case respectively. Explicitly including 

boundary cases ensures that the sets Bi and Bi^n are disjoint for all i and (/,n), i.e., these sets are either 
open or half-open sets such that no two sets share a power policy in common. This in turn simplifies the 
convex optimization as follows. 

Let P^^^ {H_) be the optimal policy maximizing the sum-rate for case i over all P_{H_) G V. The optimal 
P_^'^\H_) must satisfy the conditions for case i, i.e., P^\H) £ Bi. If the conditions are satisfied, we prove 
the optimality of P^'^\iP) using the fact that the rate function for each case is concave. On the other 
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Case 3a 



Case 36 



Case 3c 



R, 




Ri 



Ri 




Fig. 3. Rate regions Rr{P) and Rd{P) and sum-rate for cases 3a, 3&, and 3c. 




- Relay 

- Best. 



Case 1 & Case 3a 



Case 1 & Case 3b 



Case 1 & Case 3c 



Fig. 4. Rate regions Rr{P) and Rd{P) for cases (1,3a), (1,3b), and (1,3c). 



hand, when P_^^\H_) ^ Bi, it can be shown that Ri + R2 achieves its maximum outside Bi. The proof 
again follows from the fact that Ri + R2 for all cases is a concave function of P,{H_) for all P_{H) € 
V. Thus, when P^'^Hh) Bi, for every P{H_) € Bi there exists a P'(H) G Bi with a larger sum-rate. 
Combining this with the fact that the sum-rate expressions are continuous while transitioning from one 
case to another at the boundary of the open set Bi, ensures that the maximal sum-rate is achieved by 
some PJyiP) ^ Bi. Similar arguments justify maximizing the optimal policy for each case over all V. 
Due to the concavity of the rate functions, only one Py\iP) or Pj''''"'\H_) will satisfy the conditions for 
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rymin rj 




R. 



— Relay 

— Dest. 



^{2},r ~ ^{2}:d 

rymin r> 

^{l},d ~ ^{l},r 



R 



Case 2 & Case 3a 



Case 2 & Case 3b 



Case 2 & Case 3c 



Fig. 5. Rate regions Rr{P) and Rd{P) for cases (2,3a), (2,3b), and (2,3c). 



its case. The optimal P*{H) is given by this P^^{IP) or 

The optimization problem for case i or case (Z,n) is given by 



maxD^ ' or max 

P€V P€V 



S.t. E[Pk{H)]<Pk,k = l,2,r 
Pk{H) > 0, k = 1,2, r 



(23) 



where 



5(1) = 


+ -^{2},r 




5(2) = 


-R{l},r + -^{2},d 




5(^) = 


RiCj for = 


(3a, r) , (36, d) 


S(3c) __ 


= RK,r S.t. i?/c,r = 


= RK.,d- 


g{l,n) 


= 5(0 s.t. 5(0 = 


5("). 



(24) 

The optimal policy for each case is determined using Lagrange multipliers and the Karush-Kuhn-Tucker 
(KKT) conditions [32, 5.5.3]. A detailed analysis is developed in the Appendix and we summarize the 
KKT conditions and the optimal policies for all cases below. From ((24l) . the KKT conditions for each 
case X, X = i, {l,n), for all i and {l,n) is given as 



/f ) {P{h)) - z^fc In 2 < 0, with equality for P^ (h) > 0, = 1, 2, i 



(25) 
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where Uk, for all k = 1,2, r, are dual variables associated with the power constraints in (I23l ). Specializing 
the KKT conditions for each case, we obtain the optimal policies for each case as summarized below 
following which we list the conditions that the optimal policy for each case needs to satisfy. 
Case 1 : The functions f[^^ {Eih)), k = l,2,r, in ^ for case 1 are 

4" (£(ft)) = (.^i^^'-.yj;^,,) (t.m) = (l.rf),(2,r) (26) 

It is straightforward to verify that these KKT conditions simplify to 

(/^) = (i;ri2-CT)^ (fc,m) = (l,d),(2,r) (28) 

and 

Case 2 : From (|24l) . since 5^^^ can be obtained from S^^^ by interchanging the user indexes 1 and 
2, the functions {P.{!l)), and hence, the KKT conditions for this case can be obtained by replacing 
the superscript (1) by (2) and using the pairs {k,m) = {l,d), (2,r) in (I26l)-(l28l). The resulting optimal 
policies are 

pP (h) = (h) for all h, and 

Pi'^ (l^) = {ldr2 - JKZ^y {k,m) = {l,r),{2,d). ^^^^ 
Case 3a : The functions fj^^""^ {P.{!l))' ^ = li2, satisfying the KKT conditions in (l25l) are 

ft^ {Em = \hr,k\^ / (l + E \hr,k\^ Pk (h)/ e'j k = l, 2. (31) 

Since this case maximizes the multiaccess sum-rate at the relay, the optimal user policies are multiuser 
opportunistic water-filling solutions given by 

\h..i\^ ^ p(3a) _ f e e \+ p(3a)_^ 



(32) 



Ife^.iT < p(3a) ^ ^ „ p{3a) 

where without loss of generality, the users are time-duplexed even when their scaled fading states in (1321) 
are the same. While the relay power does not explicitly appear in the optimization, since this case results 
when the sum-rate is smaller than that at the destination, choosing the optimal relay policy to maximize 
the sum-rate at the destination, i.e., P^'^"-* (H) = Pr^'* (P_), suffices. 

Case 36 : The functions /^^^^ {E{!l))^ ^ = 1)2, satisfying the KKT conditions in (l25l ) can be obtained 
from (EB by replacing the subscript V by 'd' in dlD while /i^^^ (P(/i)) = (P(/i)) = /i^^ (P(/i))- 
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Thus, this case maximizes the multiaccess sum-rate at the destination and the optimal user policies are 
multiuser opportunistic water-filling solutions given by 



(33) 



p(3'^)(/,) = f^_ p(3a)^ 

Vl V2 1 \ui\n2 \hdi\ J ' 2 

p(3a)(;^)^0p(3a)^/ , _ , \ + 
Vl — U2 1 ' ^ \ 1^2 In 2 \hd,2\ J 

while the optimal relay policy is a water-filling solution P^^^"-* [H) = f^^^^ (H). 

Case 3c : The functions /^'^^^ {P.{h)), k = l,2,r, satisfying the KKT conditions in (|25l ) are given as 

/f ^ = (1 - a) /f + a/f ) , A; = 1, 2, (34) 

(P(A)) = a/(=^') (Em , k = r, (35) 
where the Lagrange multiplier a accounts for the boundary condition 

Rjc,diPiH)) = R^,r{P{H)) (36) 
and the optimal policy (H) G B^c satisfies this condition where B^c is the set of P {H_) that satisfy 



(I36b . Thus, this case maximizes the multiaccess sum-rate at both the relay and the destination. In the 
Appendix, using the KKT conditions we show that the optimal user policies are opportunistic in form 
and are given by 

ft^/^i > ^V^2 (h) = (root of ^2=0)^ (h) = 



/f V^i < ft^H Pi''^ (h) = 0,Pt^ [h] = (root of ^)|p.=o' " ^^^^ 



Analogous to cases 3a and 35, the scheduling conditions in (1371) depend on both the channel states 
and the water-filling levels i/^ at both users. However, the conditions in ( [37] ) also depend on the power 
policies, and thus, the optimal solutions are no longer water-filling solutions. In the Appendix we show 
that the optimal user policies can be computed using an iterative non-water-filling algorithm which starts 
by fixing the power policy of one user, computing that of the other, and vice-versa until the policies 
converge to the optimal policy. The iterative algorithm is computed for increasing values of a € (0, 1) 
until the optimal policy satisfies (l36l ) at the optimal a*. The proof of convergence is detailed in the 
Appendix. 

Boundary Cases (/,n) : A boundary case (/,n) results when 

5(0 ^ gin) 1 = 1^2, and n = 3a, 35, 3c. (38) 

Recall that S^^'^ and 5''^"^ are sum-rates for an inactive case /, and an active case n, respectively. Thus, 
in addition to the constraints in ( |23] ). the maximization problem for these cases includes the additional 
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constraint in (l38l ). For all except the two cases where n = 3c, the equality condition in (l23l ) is represented 
by a Lagrange multiplier a. The two cases with n = 3c have two Lagrange multipliers ai and Q2 to 
also account for the condition S^^""^ = S^^^\ 

For the different boundary cases, the functions /^^'"^ iP.{h)), k = 1,2, satisfying the KKT conditions 



in ( [25] ) are given as 

/f (Pih)) = (!-«) /i'^ + a/i"^ (P(/l)) , = 1, 2, n / 3c (39) 

/f 3^=) = as) /^'^ (P(/i)) + as/f ^ {m)) + ^ (Z(Z^)) , k = l,2, (40) 

fr'"\P{h)) = afP {P{h)) , n = 3a (41) 

(£(/i)) = a/J') (P(/i)) + (!-«) (£(/i)) , n = 36 (42) 

= ai/i') (P(/i)) + (1 - ai - a2) /i''^ , n = 3c (43) 



For ease of exposition and brevity, we summarize the KKT conditions and the optimal policies for case 
(1,3a). In the Appendix, using the KKT conditions we show that the optimal user policies pI^'^"'\h) 
are opportunistic in form and are given by 

> ^ Pi{h) = (root of F^^'^'^) |p,=o) , P2 (/i) = 

,(l,3a) r(l,3a) , , / i ^ ■\n\ . \+ ' 

^ < ^ Pi (Zl) = 0, P2 (/i) = (root of P^^^-'") |p,=o) 
where F^^''^"^ = /^^'^"^ — Ukhi2, k = 1,2. As in case 3c, the optimal policies take an opportunistic 
non-waterfilling form and in fact can be obtained by an iterative non-water-filling algorithm as described 
for case 3c. The optimal pj^'^"^ (H) = aPj^^ {W) is a water-filling solution. 

The optimal policies for all other boundary cases can be obtained similarly as detailed in the Appendix. 
In general, for all boundary cases, the optimal user policies are opportunistic non-water-filling solutions 
while that for the relay are water-filling solutions. Finally, the sum-rate maximizing policy for any case 
is the optimal policy only if it satisfies the conditions for that case. The conditions for the cases are 



Case 1 : P^^^^ < P-^^-^ and P{2},r < Pfsf,^ ^^S) 

Case 2 : R^,^^^^. < R^^^^ and P{2},d < Pfa^r ^^^^ 

Case 3a : (P/c), < (P/c)^ (47) 

Case 35 : (P^;), > (P/c)d (48) 

Case 3c : (P^)^ = (P^)^ (49) 
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Case 


(1, 3o) 


RK,r 


= R{1}4 


+ i?{2},, < {RK)d 


(50) 


Case 


(2,3a) 


RlC,r 


= R{l},r 


+ R{2},d < {RK)d 


(51) 


Case 


(1,36) : 


RK.,d 


= Ru\ d 


+ R{2},T < iRjc)r 


(52) 


Case 


(2,36) : 


RK,d 


= R{l},r 


+ R{2},d < {RK)r 


(53) 


Case 


(1,3c) : 


RK.,r 


= RlC,d = 


- R{l},d + R{2},r 


(54) 


Case 


(2,3c) : 


RK,r 


= Rk4 - 


= R{l},r + R{2},d 


(55) 



where in fading state H, are evaluated for CM \^,Pl''' (H) /ej,k = l,2, and ~ 

CM (o, Pj^^ (H) /f^ for X = i, (/, n). 

The following theorem summarizes the form of P* and presents an algorithm to compute it. 

Theorem 3: The optimal policy P_*{H_) maximizing the DF sum-rate of a two-user ergodic fading 
orthogonal MARC is obtained by computing P^^\H) and P^^'^\lP) starting with the inactive cases 1 
and 2, followed by the boundary cases (/,n), and finally the active cases 3o, 36, and 3c until for some 
case the corresponding P_^'^\H_) or P_^^''^\H_) satisfies the case conditions. The optimal P*{H_) is given 
by the optimal P^*) {H_) or p(^'") (^H_) that satisfies its case conditions and falls into one of the following 
three categories: 

Inactive Cases: The optimal policy for the two users is such that one user water-fills over its link to 
the relay while the other water-fills over its link to the destination. The optimal relay policy P^{iP) is 
water-filling over its direct link to the destination. 

Cases (3a, 36, 3c): The optimal user policy Pj^{H_), for all € /C, is opportunistic water-filling over its 
link to the relay for case 3a and to the destination for case 36. For case 3c, P^{H_), for all k £ IC, takes 
an opportunistic non-waterfilling form and depends on the channel gains of user k at both receivers. The 
optimal relay policy P* {H_) is water-filling over its direct link to the destination. 

Boundary Cases: The optimal user policy P^{H_), for all k € )C, takes an opportunistic non-water- 
filling form. The optimal relay policy P* {H_) is water-filling over its direct link to the destination. 

Proof: The closed form expressions for the optimal policies for each case are developed in the 
Appendix. The need for an order in evaluating P_*{H_) is due to the following reasons. Since every case 
results from an intersection of two polymatroids, the conditions in (I47l)-(l49l) hold for all cases. Thus, 
all feasible power policies satisfy one of these three conditions as a result of which these conditions 
do not allow a clear distinction between the cases. In contrast, the conditions for cases 1 and 2 in (l45l) 
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and (l46l ). respectively, are mutually exclusive. For the boundary cases, since every boundary case (/,n) 
results from the intersection of an active case n = 3a, 36, 3c, with an inactive case Z = 1, 2, and is itself 
an active case, one of its conditions corresponds to the condition for case n, n = 3a, 36,3c in (|47])-(|49]). 
Additionally, in the Appendix we show that the boundary condition S^^^ = 5"^"^ implies that only one 
of the two inequality conditions of case I holds strictly while the other simplifies to an equality. An 
immediate implication of these two conditions is that the boundary cases are mutually exclusive and 
the set of power policies satisfying them are also disjoint from those satisfying cases 1 and 2. Thus, to 
determine the optimal P_*{H_), one can start with any one of the mutually exclusive inactive and boundary 
cases. If the optimal policy for any one of these cases satisfies its case conditions, then, P_*{H_) is given 
by that policy. However, if all these cases are eliminated, i.e., none of their optimal policies satisfy the 
appropriate case conditions, the optimal policies for remaining three cases 3a, 36, and 3c can be computed 
one at a time. From (|47])-(|49l). cases 3a, 36, and 3c are mutually exclusive, i.e., their feasible power sets 
are disjoint, and thus, the optimal policy, satisfies the conditions for only one of three active cases. ■ 

Remark 5: The conditions for cases 3a, 36, and 3c can also be redefined to include the negation of all 
the conditions for the other cases. This in turn eliminates the need for an order in computing the optimal 
policy; however, the number of conditions that need to be checked to verify if the optimal policy satisfies 
the conditions for cases 3a or 36 or 3c remain unchanged relative to the algorithm in Theorem [3] 

We now discuss in detail the optimal power policies at the sources and the relay for the different cases. 

Optimal Relay Policy: In the orthogonal model we consider, the relay transmits directly to the desti- 
nation on a channel orthogonal to the source transmissions. Thus, the relay to destination link can be 
viewed as a fading point to point link. In fact, in all cases the optimal relay policy involves water-filling 
over the fading states analogous to a fading point to point link (see [33]). However, the exact solution, 
including scale factors, depends on the case considered. 

The optimal cooperation strategy at the relay also depends on the case studied. For instance, consider 
case 1 where users 1 and 2 achieve significantly larger rates (relative to the other receiver) at the relay 
and destination, respectively. Thus, the sum-rate is the sum of the rates achieved over the bottle-neck 
links from user 1 to the destination and from user 2 to the relay; i.e., it is the sum of the single-user rate 
user 1 achieves at the destination and the rate user 2 achieves at the relay. The single-user rate achieved 
by user 1 at the destination requires the relay to completely cooperate with user 1, i.e., the relay uses 
its power Pr {H_) to forward only the message from user 1 in every fading state in which it transmits. 
As shown in the intersection for case 1 in Fig. |2l this is due to the fact that since user 1 achieves a 
significantly larger rate at the relay than does user 2, the sum-rate is maximized when the relay allocates 
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its resources entirely to cooperating with user 1. Finally, for case 2, the relay cooperates entirely with 
user 2. 

For the active cases, 3a and 36, the sum-rate may be achieved by an infinite number of feasible points 
on one or both of the sum-rate planes; the optimal cooperative strategy at the relay will differ for each 
such point. Thus, for a corner point the relay transmits a message from only one of the users while for 
all non-comer points the relay transmits both messages. 

For the boundary cases including case 3c, the requirement of an equality (boundary) condition results 
in the introduction of an additional parameter. Thus, for case 3c, the parameter a is introduced to satisfy 
the equality constraint on the sum-rates achieved at the relay and destination. Similarly for cases (1,3a), 
(2,3a), (1,36), and (2,36), the parameter a is chosen to ensure that the optimal power policies at the 
users and relay satisfy the equality constraint for that case. Finally, for cases (1,3c) and (2,3c), the 
requirement of satisfying two boundary conditions requires the introduction of the two parameters, a\ 
and 02, one for each condition. 

Optimal User Policies: As with the relay, the optimal policies for the two users depend on the case 
considered. For cases 1 and 2, the optimal policies are water-filling solutions, i.e., each user allocates 
power optimally as if it were transmitting to only that receiver at which it achieves a lower rate. In 
fact, the conditions for case 1 in ( |45l) suggest a network geometry in which source 1 and the relay are 
physically proximal enough to form a cluster and source 2 and the destination form another cluster; 
and vice-versa for case 2. This clustering and the resulting water-filling over the bottle-neck links is the 
reason why the relay forwards the message of only that user physically proximal to it, namely, only user 
1 and only user 2 for cases 1 and 2, respectively. 

For case 3a, the optimal policies at the two users maximize the sum-rate achieved at the relay (the 
smaller of sum-rates achieved at the two receivers). These policies are the same as those achieving the 
sum-capacity of a two-user multiple-access channel with the relay as the intended receiver (see [27], 
[28]). Thus, the optimal policy for each user involves water-filling over its fading states to the relay. The 
solution also exploits the multiuser diversity to opportunistically schedule the users in each use of the 
channel. 

Analogously, the optimal policies for case 36 require multiuser water-filling over the user links to the 
destination. For both cases, if the channel gains have a joint fading distribution with a continuous density, 
the sum-rate maximization simplifies to scheduling only one user in each fading instantiation (parallel 
channel). Thus, the users time-share their channel use and the maximum sum-rate is achieved by a unique 
point on the boundary of the rate region (see [28, III.D]). 
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The optimal policies for case 3c require the users to allocate power such that the sum-rates achieved 
at both the relay and the destination are the same. This constraint has the effect that it preserves the 
opportunistic scheduling since the sum-rate involves the multiaccess sum-rate bounds at both receivers. 
However, the solutions are no longer waterfilling due to the fact that the equality (boundary) condition 
results in the function fjf'^^ being a weighted sum of the functions fj^^""^ and fj^^^^ for cases 3a and 36, 
respectively. 

Finally, the requirement of satisfying one or more boundary conditions also affects the nature of the 
optimal policies for all the (/,n) cases. Thus, for these boundary cases, since the sum-rate planes are 
active, i.e., the functions fj:''^^ involve the multiple-access sum-rate bounds, the optimal power policies 
result in an opportunistic scheduling of the users. However, as with case 3c, here too the optimal policies 
are no longer water-filling since the boundary conditions result in the functions fj!'"^ being weighted 
sums of the functions for cases I and n. 

Remark 6: The case conditions in (|45])-(|55]) require averaging over the channel states; thus, the case 
that maximizes the sum-rate depends on the average power constraints and the channel statistics (including 
network topology). 

Remark 7: The optimal policy for each source for cases 1, 2, 3a, and 36 depends on the channel 
gains at only one of the receivers. However, the optimal policy for the boundary cases, including case 
3c, depends on the instantaneous channel states at both receivers. Furthermore, all the cases exploiting 
the multiuser diversity require a centralized protocol to coordinate the opportunistic scheduling of users. 

V. Two-User DF Rate Region: Optimal Power Policies 
In Theorem |2j the DF rate region TZdf is shown to be a union of the intersections of the regions 
T^r{P.iIL)) and TZdiE.ilL)) achieved at the relay and destination, respectively, where the union is over 
all P(H) G V. Furthermore, since TZ^p is convex, each point on the boundary of TZdf is obtained by 
maximizing the weighted sum fiiRi + IJ.2R2 over all P_{H_) G V, and for all fii > 0, fi2 > 0. In fact, 
for every {111,112), the rate tuple (i?i,i?2) that maximizes the weighted sum ^iiRi + ^2-^2 results from 
an intersection of two rate polymatroids. 

Thus, analogously to the sum-rate analysis for fJ-i = fi2 = 1, for arbitrary (/Ui,/X2), ^i-Ri + ^2^2> is 
maximized by either one of two inactive cases, or by one of nine active cases of which six are boundary 
cases. To find the rate tuple maximizing fiiRi + fi2R2, we use the classic result in linear programming 
that the maximum value of a linear function constrained over a feasible bounded polyhedron is achieved 
at a vertex of the polyhedron [32, Chapter 1.2.2]. Thus, for any P_{iP), the (i?i, i?2)-tuple maximizing 
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fiiRi + /U2-R2 is given by a vertex of a TZr {P.{H_)) n TZd {E.{IL)) polyhedron at which fiiRi + 112R2 is 
a tangent. Recall that for the two inactive cases, the polymatroid intersections result in rectangles, and 
thus, there is a unique vertex maximizing /ii-Ri + /i2-R2- The intersection is also a rectangle for the six 
boundary cases since these active cases are such that only one point on the sum-rate plane is included in 
the region of intersection. On the other hand, for cases 3a, 36, and 3c, the intersection of i^-dimensional 
polymatroids results in a K-dimensional polyhedron. Thus, for these cases, when /^i < /X2, ^iiRi + ^2^2 
is maximized by that vertex where user 1 is decoded before user 2, i.e., at the vertex where user 2 
achieves its maximal single-user rate. 

For simplicity, we present the results for cases 1, 3a, and (1, 3a) . The results for the other cases follow 
naturally from discussions for these three cases. Without loss of generality, we let fii < the analysis 
for /xi > /X2 follows in an analogous manner 

Case 1: From Fig. |2j the weighted sum fiiRi + ^2^2 for this case is given by 

f^lR{l},d + f^2R{2},r- (56) 

Since fii and fi2 are independent of the transmit powers, the optimization problem is the same as that 
for the sum-rate case. Thus, at the maximal rate, users 1 and 2 waterfill on their bottleneck links to the 
destination and relay, respectively. 

The analysis for case 2 is the same as that for case 1 except now user 1 and 2 waterfill on their 
bottleneck links to the relay and destination, respectively. 

Case 3a: The weighted sum fiiRi + /U2-R2 is maximized by the vertex with coordinates 

Rl = RK,r - R2 (57) 

R2 = min {R{2},r,R{2},d) ■ (58) 
The maximization of + /X2-R2 thus simplifies to 

max {fiiRic,r + (/U2 - ^1) min {R{2},r, R{2},d)) (59) 

where the optimal P^^"-\h_) satisfies the conditions in (|47] ) for this case. As in the appendix, there are 
three possible disjoint solutions to the max-min optimization in ( [59l ) resulting from either R{2}.r being 
smaller, larger, or equal to R{2},d- We discuss the optimal policies for each of these sub-cases separately 
below. 

1) ^{2},r < R{2},d- For this sub-case, the vertex of interest in (l57l) and (l58l ) is achieved by the 
MAC bounds at the same receiver Thus, the maximization for these cases simplifies to that 
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developed in for an ergodic fading MAC in [28]. The optimal power policies involve water-filling 
and opportunistic scheduling of the users and water-filling at the relay over its direct link to the 
destination. 

2) R{2},r > R{2},d- The maximization here simplifies to 

max {fiiRjc,r + (M2 - /^i) R{2},d) ■ (60) 

As with the Lagrangian expressions for the boundary cases in the Appendix, here too, the weighted 
sum of rates in (l60l ) is an appropriately weighted mixture of sum and single-user rates achieved at 
the relay and destination, respectively. Thus, analogously to the boundary cases, one can verify in a 
straightforward manner that the optimal policies maximizing (l60l) at both users are non-waterfilling 
solutions with opportunistic scheduling based on relative fading states while that at the relay requires 
waterfilling over its direct link to the destination. Note that the optimal policies at both the users 
and the relay depend on the values chosen for /xi and fi2. 

3) R{2}.r = R{2},d- Subject to average power and positivity constraints, the maximization here 
simplifies to 

maxpgp {fllRjC,r + (/^2 - m)-R{2},r) 

(61) 

S.t. i?{2},r = R{2},d- 

The maximization in (|6TI ) subject to the equality constraint results in a Lagrangian with a weighted 
mixture of single-user rates achieved at the relay and destination. Thus, the optimal user and relay 
policies have a form similar to that discussed in the previous sub-case in which the sum-rate and 
single-user rate are achieved at different receivers. 
Remark 8: The three sub-cases for case 3a studied above are differentiated by additional constraints 
relating the single-user rates at the relay and the destination. This in turn implies that the region B^a 
will be divided into three mutually exclusive subsets, where the condition for each sub-case is satisfied 
in only one of the subsets. 

Boundary case (1,3a): Recall that a boundary case {l,n) results from satisfying the conditions for 
the active case n and satisfying the conditions for the inactive case I as a mixture of equalities and 
inequalities. The resulting rate region belongs to the set of active cases but has one unique sum-rate 
point such that the intersection of pentagons results in a rectangle (see Figs. |4] and |5]). Thus, the weighted 
optimization + /U2-R2 for case (1,3a) simpUfies to 

fJ'lR{l},d + ^J'2R{2},r 

(62) 

s.t. (i?/c),. = R{i}^d + R{2},r- 
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Note that the constraint in (l62l ) is the same as that for the boundary case (1,3a) in (l50l ). Thus, the 
constrained maximization problem in (l62l ) is analogous to the sum-rate maximization for the boundary 
cases and admits a similar non-waterfilling opportunistic solution for the user power policies and a 
waterfilling solution at the relay. 

Remark 9: The discussion here for cases 3a and (1,3a) also applies to the other active (including 
boundary) cases. In each such case, the optimal policies depend on all the Lagrange dual variables, with 
each variable reflecting a specific constraint. 



A. K-user Sum-Rate Analysis 

We use Lemma [T] to extend the analysis in the previous sections to the iT-user case. Recall that TZdf 
is given by a union of the intersection of polymatroids, where the union is over all power policies. From 
Lemma [T] we have that the maximal K-user sum-rate tuple is achieved by an intersection that either 
belongs to active set or to the inactive set. We write / = 1, 2, . . . , 2^ — 2, to index the 2^ — 2 non-empty 



VL /C-USER Generalization 




(63) 




= i?yc,i for (i, j) = (3a, r) , (36, d) (64) 
S^^"^) = R,c,r s.t. R,c,r = RK,d- (65) 
Finally, the sum-rate for the boundary cases totaling 3 (2^ — 2) and enumerated as cases {l,n), 



/ = 1, 2, . . . , 2^ - 2, n = 3a, 36, 3c, are 

case {I, 3a) : SC'^'^) = RK,r 




(67) 



s.t. Rjc4 = Rs,r + RK.\s,d for case / 

case (/,3c) : ^(''3^) = RK,r 

s.t. Rjcr = Rk.4 = Rs,r + RK:\s,d for case / 



February 10, 2009 



DRAFT 



25 



where the subset S is chosen to correspond to the appropriate case I. 

Remark 10: The constraint for case / in ( [63] ) can be obtained directly from the requirement that the 
K-user sum-rate constraints at the two receivers are larger than that for case I (see ([8]l). 

The K-user sum-rate optimization problem for cases i and {l,n) can be written as 



max ' or max 

Per Pev 

s.t. E[Pk{H)]<Pk,k = l,2,r 

Pk{H)>0, k = l,2,r. I ^^g^ 

An inactive case / results when the conditions for that case in (1631 ) are satisfied. A boundary case results 
when the conditions for one of the cases in (l66l)-(l68]) is satisfied for the appropriate (/,n) case. Finally, 
case 3a or 36 or 3c results when the conditions for neither the inactive nor the boundary cases are 
satisfied. 

As in Section |lVl the optimization for each case involves writing the Lagrangian and the KKT 
conditions. The optimal policy P_*{H_) satisfies the conditions for only one of the cases. For brevity, 
we summarize the details below. 

• Inactive cases: The Lagrangian for these cases involves a sum of the DF bounds at the relay in ([T8] ) 
for users in S, for each non-empty S, and the DF bounds at the destination in ( [T9l) for the remaining 
users in }C\S such that for i = 1, 2, . . . , 2^ — 2, 



where u^, for all k, are the dual variables associated with the power constraints in and > 
are the dual variables associated with the positivity constraints (P^ > 0) on P^. Writing the KKT 
conditions, it is straightforward to verify that the optimal policy for a user in 5 is a function of the 
channel gains at the relay while that for a user in }C\S is a function of the channel gains only at 
the destination. In fact, when S or }C\S are singleton sets, the optimal policy for the user in S or 
JC\S is simply water-filling over its bottle-neck link to either the relay (if S) or the destination (if 
}C\S). More generally, when S or JC\S are not singleton sets, the optimal policy is an opportunistic 
water-filling solution. Finally, the relay's policy is water-filling over its direct link to the destination. 
Cases 3a, 36, and 3c: The Lagrangian for these three cases is given by 

£W= S^'^ -ZkeT^k^[Pk{H)-Pk]+EkeT^kPk{H) i = 3a, 36, 3c (71) 
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where 



E 
E 



i = 3a 

i = 36 (72) 



oc [Ekeic\Hr,k\ Pkim/e 

a5(3«) + (1 - a) ^ = 3c 

where the dual variable a is associated with the boundary condition = S^^''^ for case 3c. From 
(TtTI) and (1721 ). for case 3a, since the dominant bounds are the MAC bounds at the relay, the optimal 
user policies involve opportunistic water-filling over their links to the relay. The optimal policies 
take a similar form for case 36, except now since the dominant bounds are the MAC bounds at the 
destination, each user opportunistically waterfills over its link to the destination. Finally, for case 
3c, the KKT conditions satisfied by Pk{H_), for all k, are 



a\hr,k\^ I {l-a)\ha_k 



YTuEn + ^^^^ 'i2P. — < with equality if Pk{h) > 0. (73) 



Thus, as in the two-user analysis for case 3c in the Appendix, the optimal user policies are no longer 
water-filling but involve opportunistic scheduling of the users to exploit the multiuser diversity. In 
fact, the optimal policy for each user depend on its channel gains to both the relay and the destination 
and can be computed using the iterative algorithm detailed in the Appendix. Finally, for all cases, 
the optimal relay policy is a water-filling solution. 
Boundary cases {l,n): The Lagrangian for these cases is given by 

I = l,2,n = 3a,3b (74) 

jril,3c) ^ ^^gil) ^ ^^giSa) + (l _ „^ _ ^(Sb) _ ^^^^ [p^(^) _ p^] 

+ EkeT>^kPkiH), 1 = 1,2 (75) 

where a is the dual variable associated with the boundary condition 5^'^ = 5*^"^ for n 7^ 3c and ai 
and a2 are dual variables associated with the boundary conditions S*-'^ = S*-^"-* and S^^"'^ = S^^''\ 
Here again the optimal solution for each is no longer water-filling and depends on the channel 
gains to both the relay and the destination. Furthermore, as with case 3c, here too the optimal user 
policies exploit the multiuser diversity to opportunistically schedule the user transmissions. Finally, 
the optimal relay policy for all boundary cases is a water-filling solution over its direct link to the 
destination. 
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Theorem 4: The optimal power policy P* {H) that maximizes the DF sum-rate of an K-user ergodic 
fading orthogonal Gaussian MARC is obtained by computing P^*^ {W) and pC'") (^H) starting with the 
inactive cases 1, 2, . . . , 2^ — 2, followed by the boundary cases (/,n), and finally the active cases 3a, 
36, and 3c until for some case the corresponding (H) or P(^'")(^) satisfies the case conditions. The 
optimal P* {H) is given by the optimal P^*-* (H) or p(''") (^H) that satisfies its case conditions and falls 
into one of the following three categories: 

Inactive Cases: The optimal user poUcy P^{H), for all A; G /C, is multi-user water-filUng over its bottle- 
neck (rate limiting) link to the relay or the destination. The optimal relay policy P* {H_) is water-filling 
over its direct link to the destination. 

Active Cases (3a, 36, 3c): The optimal user policy P*„(H), for all /c € /C, is opportunistic water-filling 
over its Unk to the relay for case 3a and to the destination for case 36. For case 3c, P^{H_), for all 
k e IC, takes an opportunistic non-waterfilling form. The optimal relay policy P* is water-filling over 
the relay-destination link. 

Boundary Cases: The optimal user policy P^{l£), for all k e JC, takes an opportunistic non- water- 
filling form. The optimal relay poUcy P* {H_) is water-filling over its direct link to the destination. 

Based on the optimal DF policies, one can conclude that the topology of the network affects the 
form of the solution with the classic multiuser opportunistic waterfilling solutions applicable only for the 
sources-relay or the relay-destination clustered models. For all other partially clustered or non-clustered 
networks, the solutions are a combination of single- and multi-user water-filling and non-waterfilling but 
opportunistic solutions. 

B. K-user Rate Region 

Analogous to the two-user analysis, one can also generalize the sum-rate analysis above to develop 
the optimal policies for all points on the boundary of the K-user DF rate region. For brevity, we outUne 
the approach below. 

We start with the observation that the DF rate region, TZdf, is convex, and thus, every point on the 
boundary of TZdf is obtained by maximizing the weighted sum Y2keKl-^kRk, ^fc > for all k. As noted 
earlier, each point on the boundary of TZdf obtained by an intersection of two polymatroids for some 
P.{H_). Thus, analogously to the sum-rate analysis for fi^ = 1 for all k, for arbitrary (/ii,/i2, ■ ■ ■ jfJ-x), 
^j^.g^/ifcPjfc, is maximized by either by an inactive or an active case. 

Since the maximum value of YlkeK/^^Pk over a feasible bounded polyhedron is achieved at a vertex of 
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the polyhedron, for any P,{H_), the R2, . . . , i?A-)-tuple maximizing J^k^Kl^^Rk is given by a vertex 
of a TZr {P.{H.)) l~l (.E.{IL)) polyhedron at which '}ZkeK.^^^^k is a tangent. For the 2^ — 2 inactive 
cases, the polymatroid intersections are polytopes with constraints on the multiaccess rates of all users 
in S and K\S at the relay and destination, respectively. Since bounds on the multiaccess rates of / users 
result in a polymatroid with l\ vertices, the intersection of the two orthogonal sum-rate planes will result 
in a polytope with (|5|!) (|/C\5|!) vertices of which an appropriate vertex will maximize J2kelC^^kRk^ 



vertices since these active cases are such that only one point on the sum-rate plane is included in the 
region of intersection. Finally, for cases 3a, 36, and 3c, the intersection of iC-dimensional polymatroids 
results in a K-dimensional polyhedron. 

In general, the intersection of two polymatroids is not a polymatroid, and thus, unlike polymatroids, 
greedy algorithms do not maximize the weighted sum of rates. This in turn implies that closed form 
expressions are not in general possible and determining the optimal power policies requires convex 
programming techniques. However, for specific clustered geometries, we present closed form results. 

We write fi to denote a vector of weights with entries /z^, for all k. Let vr be a permutation corre- 
sponding to a decreasing order of the entries of jj, such that tt {k) is the /c*'^ entry of vr and ir {j : k) = 
{vr (j) , vr (j + 1) , . . . , vr (k)}. Thus, X^feeC^ (^) -^7r(fc) is maximized by a vertex whose rate tuple i?2, 
is such that > RTr(2) > • • • > Rtt(k)> i-e-> the decoding order at the vertex is the reverse of the 

order of entries of fi. 

For simplicity, as with the two-user analysis, we summarize the results for cases /, 3a, and (/,3a), 
for all / G {1, 2, . . . , 2^ - 2}. The results for the other cases follow naturally from discussions for these 
cases. 

Case I : This case results when the sum-rate plane at the relay for users in 5 C /C intersects the sum- 
rate plane at the destination for the complementary users in }C\S. For a permutation vr with decreasing 
order of the entries in fi, let vr^ be the decreasing order of the entries of /j, for the users in A C IC . The 
weighted rate-sum can be expanded as 




boundary cases are also characterized by an intersection with (|5|!) (|/C\5|!) 



(76) 



where (1761 ) is maximized by choosing the rates i?7r(fc) for all k as 




(77) 




(78) 
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and 



R 



R 



'Tric\s{k) 



E 
E 



C 















2 



k = 2,3,..,\IC\S\ 



(79) 
(80) 



Thus, the users in S and /C\5 are decoded in the increasing order of their weights at the relay and 
destination respectively. The optimal power and rate allocation for the users in S and }C\S are the 
multiuser opportunistic water-filling solutions at the relay and destination, respectively, and can be 
computed using a utility function approach developed in [28, II.C]. 

Case 3a: The poly tope resulting from the intersection of two polymatroids is defined by the constraints 



Rs < min|E 
and Rk<E 



kes 



C 



Pk 



Pk 



,E 



es \^d,k 



Pi 



for all 5 C /C (81) 
(82) 



However, since the polytope given by (I8TI) and (182] ) above is in general not a polymatroid, greedy 
algorithms cannot be used to maximize the weighted sums and thus developing closed form solutions 
for this case is not possible in general. However, the optimal policies maximizing the weighted sum of 
rates can be computed in strongly polynomial timen [26, Theorem 47.4]. 

Remark 11: For the special case where the bounds at the relay are smaller than the bounds at the 
destination for all S, i.e., IZr C TZd, the optimal user policies are multiuser water-filling solutions 
developed in [28, II.C] with the relay as the receiver. Note that this condition implies that all possible 
subset of users achieve better rates at the destination than at the relay. This can happen when either 
all users are clustered closer to the destination or when the relay has a relatively high SNR link to the 
destination sufficient enough to achieve rate gains for all users at the destination. 

Remark 12: Similarly, for case 36, for the special case in which TZ^ C TZr, the optimal user policies 
are multiuser water-filling solutions with the destination as the receiver. In the following section we show 
that DF achieves the capacity region when case 36 holds for all points on the boundary of the outer 
bound rate region. In fact, this condition implies that all possible subset of users achieve better rates at 
the relay than they do at the destination which in turn suggests a geometry where all subsets of users are 

^An algorithm is said to run in strongly polynomial time when the algorithm run time is independent of the numerical data 
size and is dependent only on the inherent dimensions of the problem. In contrast polynomial time algorithms are characterized 
by run times that are polynomial not in the size of the input but the numerical value of the input which may be exponentially 
large. 
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clustered closer to the relay than to the destination. The optimal relay policy in all cases is a waterfilling 
solution over its link to the destination. 

Boundary case (Z,3a): Recall that a boundary case (/,n) results when the K-mtx sum-rate for the 
active case n is equal to that for the inactive case /. The resulting region of intersection, analogous to 
the inactive cases, is a polytope with (5!) (/C\5!) vertices. The weighted optimization "^^(z^qT^ (k) R-j^i^k) 
for case (/,3a) simplifies to 

(oj) 

s.t. {Ric)r = T^keK^ (^) ^7r(fe) = T^kes^TTsik) + T.kes\!C^7T{k) 
where -R7rs(fe) RnK:\s(k) given by dTTt-dSOl). Here again, given the complexity of the optimization, 
closed form solutions are difficult to obtain. However, as before, one can compute the optimal policies 
and the rate tuple maximizing ([83] ) in polynomial time using combinatorial methods. 

Remark 13: The discussion here for cases 3a and (1,3a) also applies to the other active (including 
boundary) cases. In each such case, the optimal policies depend on all the Lagrange dual variables, with 
each variable reflecting a specific constraint. 

vn. Outer Bounds 

An outer bound on the capacity region Cm arc of a i('-user full-duplex MARC is presented in [10] (see 
also [29, Th. 1]) using cut-set bounds as applied to the case of independent sources and we summarize 
it below. 

Proposition 4 ( [29, Th. 1]): The capacity region Cmarc is contained in the union of the set of rate 
tuples {Ri, R2, ■ ■ ■ , Rk) that satisfy, for all S Q K, 

Rs<m.m{I{Xs;Yr,Yd\Xs^,Xr,U)J{Xs,Xr;Yd\Xs^,U)} (84) 

where the union is over all distributions that factor as 

p{u) ■ (^Ilfcli^^^^l^)) -pixrlxfcu) ■p{yr,yd\xK.,Xr). (85) 
Remark 14: The time-sharing random variable U ensures that the region in ([84] ) is convex. One 
can apply Caratheodory's theorem [34] to this i^-dimensional convex region to bound the cardinality of 

U as \U\ <K + l. 

Proposition 5: For the orthogonal MARC the cutset bounds in ([84] ) specialize as 

Rs< min{0/(X5;i;yrf,i|X5e,C/),^/(X5;yd,ilX5.,?7)+^/(X,,;yrf,2|C/)} foraU5C/C (86) 
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where the union is taken over all distributions that factor as 



p{u) 



(87) 



nfc=i^(^'^'l")j 'Piyryd\xK.) + 9 ■p{Xr\u) ■ p{yd\Xr 

Remark 15: The above bounds can also be obtained by using a mode variable Mr to denote the 
half-duplex listen and transmit states at the relay such that Mr is in the listen and transmit states with 
probabilities 6 and 1 — 6, respectively (see [35]). The instantaneous relay mode is assumed known at all 
nodes, such that (l86l ) results from conditioning the bounds in (l84l ) on Mr, and (l87l ) from replacing Xr 
with {Xr,Mr) in dSSl ) and expanding the resulting joint distribution. 

Remark 16: The joint distribution for the cutset bounds in (|87]) is the same as that for DF in ([T2l) . 
This is in contrast to the full-duplex MARC where in general, the two distributions (and bounds) are not 
the same. 

Theorem 5: For a degraded orthogonal discrete memory less MARC where Xg—Yr—Yd form a Markov 
chain, DF achieves the capacity region of a degraded orthogonal MARC. 

Proof: The proof follows directly from applying the Markov property Xg — Yr — Yd to the cutset 
bounds in ([86l) and comparing the resulting bounds with those for DF in ([TT]) . Note that for the full- 
duplex degraded MARC, the inner and outer bounds are not the same in general. In fact, for the degraded 
Gaussian (full-duplex) MARC, it has been recently shown in [29] that DF achieves the sum-capacity when 
the intersection of the two polymatroids at the relay and destination belongs to the set of active cases. ■ 

For an orthogonal Gaussian MARC with fixed and 9, using a conditional entropy theorem, one 
can show that Gaussian signals maximize the bounds in (l86l ). Thus, substituting ~ CM {0, Pk/9), 
k = l,2, and X,. ~ CM (O, Pr/9) in ([86]l, we have 



Rs ^ mill 9 log 



I + ^GkPk/* 

fcG5 



,9c(^\Hd,k\^ Pk/9^ +9c{\Hd,rf Pr/9^^ (88) 



where 



Hr,k Hdj 



T 



^r,k ^d,k 



(89) 



and H*s^ is the complex conjugate of H(.)- Using the fact that the ergodic channel is a collection of 
parallel non-fading channels, the capacity region of an ergodic fading orthogonal Gaussian MARC is 
given by the following theorem. 

Theorem 6: The capacity region Cq-marc of an ergodic fading orthogonal Gaussian MARC is 
contained in 

noB = U {^1 (P) n 7^2 (P)} (90) 
Per 
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where, for all 5 C /C, we have 
and 

7^2 {P)=\{Ri,R2):Rs<^ 



eiog 



kes 



(91) 



(92) 



ec \^ \Hd,kf Pk{H)/ej +dc (\Hd,r\^ Prim/e'^ 

Remark 17: Comparing outer bounds in ^2) with the DF bounds in (fT9l) . we see that the bounds 
at the destination are the same in both cases. However, unlike the DF bound at only the relay in ([TSl l. 
the cutset bounds in (I9l1 ) is a SIMO bound with single-antenna transmitters and both the relay and the 
destination acting as a multi-antenna receiver. 

The expressions in (|9l1 ) and ( |92l ) are concave functions of Pk{H_), for all k, and thus, the region TZqb 
is convex. Thus, as in Theorem |2l the region TZqb in (l90l ) is a union of the intersections of the regions 
T^i{P{M.)) and 'R-2{P.{H)) , where the union is taken over all P{H_) G V and each point on the boundary 
of TZdf is obtained by maximizing the weighted sum fiiRi + ^2R2 over all P(H) G V, and for all 
fii > 0, fi2 > 0. In [29], it is shown that the rate polytopes satisfying ([84l ) are polymatroids. Since, the 
polytopes in (|9T] ) and (l92l ) are obtained from (l84l) for the special case of orthogonal signaling, one can 
verify in a straightforward manner using Definition [T] that these are polymatroids as well. 



A. Optimal Sum-rate Policies and Sum-capacity 

Since TZob is obtained completely as a union of the intersection of polymatroids, one for each choice 
of power policy. Lemma [T] can be applied to explicitly characterize the outer bounds on the sum-rate. 
Thus, the maximum sum-rate tuple is achieved by an intersection that belongs to either the active set or 
to the inactive set. Let / = 1, 2, . . . , 2^ — 2, index the 2^ — 2 non-empty subsets of K.. For a K-user 
MARC, there are (2^ — 2) possible intersections of the inactive kind with sum-rate given by 

easel : JW = /?5,i + i?x;\5,2 / = 1, 2, . . . , 2^ - 2 ^^^^ 

s.t. Rs,i < R^^ and R,c\s,2 < R^^^^i 
where Rj^j and -^^4'^ are as defined in Section JII] and for j = 1, 2, are given by the bounds in ( |9T] ) and 
(l92l ). respectively. The sum-rates J^'^\ i = 3a, 36, 3c, are 



J» = i?cj for = (3a, 1), (36, 2) 
= Rk,i s.t. Rk,i = Ric,2. 



(94) 
(95) 
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Finally, the sum-rate J^''"), for the 3 (2^ — 2) boundary cases, enumerated as cases {l,n), I 
2, n = 3a, 36, 3c, are 

case {I, 3a) : J^^'^'^) = 
s.t. Rfc,! = Rs,i + Rk.\s,2 for case I 
case {I, 3b) : J^'^^) = R,c,2 
s.t. Rk.,2 = Rs,i + Rk\s,2 for case / 
case (/,3c) : J^'^^) = R/c,! 
s.t. Rjc,i = Rk.,2 = Rs,i + Ric\s,2 for case / 
where the subset S is chosen to correspond to the appropriate case I. 
The ivT-user sum-rate optimization problem for case i and case (/, n) is 

max J^^ or max 

PeV Per 

s.t. E[Pk{H)]<Pk,k = l,2,r 
Pk{H)>0, k = l,2,r. 

An inactive case / results when the conditions for that case in ( |93] ) are satisfied. A boundary case results 
when one of the conditions in (|96l)-(l98]) is satisfied for the appropriate {I, n) case. Finally, case 3a or 36 
or 3c results when the conditions for neither the inactive nor the boundary cases are satisfied. 

As in Section |IVl the optimization for each case involves writing the Lagrangian and the KKT 
conditions. The optimal policy p!'"^\H_) satisfies the conditions for only one of the cases. For brevity 
and to avoid repetition, we summarize the details below. 

• Inactive cases: The Lagrangian for these cases involves a sum of the MIMO cutset bounds in (|9l1 ) 
for users in S, for some S, and the cutset bounds at the destination in (|92l ) for the remaining users in 
}C\S. Thus, the optimal policy for a user in 5 is a function of the channel gains at both the relay and 
destination while that for a user in 1C\S is a function of the channel gains only at the destination. For 
K > 2, using the results in [36, Theorem 1] for ergodic fading SIMO-MAC channels, the policies 
for the users in S are water-filling and allow at most = 4 users to transmit simultaneously, where 
/ is the number of antennas at the receiver. Furthermore, the optimal user policies can be obtained 
using an iterative water-filling approach [37]. On the other hand, the SISO-MAC bounds for the 
users in }C\S result in a multiuser opportunistic water-filling solution. Finally, the relay's policy is 
water-filling over its direct link to the destination. 



(96) 
(97) 
(98) 
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• Cases 3a, 36, and 3c: For case 3a, the dominant bounds are the SIMO cut-set bounds, and thus, 
as discussed for the inactive cases, the optimal poUcy is water-fiUing for each user such that a 
maximum of 4 users can transmit simultaneously. On the other hand for case 36, the dominant 
bounds are the cooperative bounds at the destination and the optimal policy for each user is an 
opportunistic water-filling solution. Finally, for case 3c, as one would expect, the optimal policies 
are no longer water-filling. In all cases, the optimal relay policy is a water-filling solution. 

• Boundary cases {l,n): The Lagrangian for these cases is a weighted sum of the sum-rates for one 
of cases 3a, 36, or 3c and one of the inactive cases. Here again the optimal solution for each is no 
longer water-filling and depends on the channel gains to both the relay and the destination. As with 
the other cases, here too the optimal relay policy for all boundary cases is a water-filling solution. 

Comparing these optimal policies with that for DF, we have the following capacity theorem. 

Theorem 7: The sum-capacity of a K-u^er ergodic fading orthogonal Gaussian MARC is achieved by 
DF when the optimal poUcy p(°'') [H) for the cutset bounds satisfies the conditions for case 36 and for 
no other case. 

Proof: The proof follows from comparing the expressions J^'^ for all cases in (l24l) and (|93])-(l98l) 
for the inner and outer bounds, respectively. For all cases where the SIMO cut-set bound dominates 
the sum-rate, the cutset bounds do not match the DF bounds. Thus, when the optimal policy (^) 
satisfies the conditions for case 36, where the sum-rate bounds at the destination dominate, DF achieves 
capacity. ■ 
Remark 18: Recall that case 36 corresponds to a clustered geometry in which the relay is clustered 
with all sources such that the cooperative multiaccess link from the sources and the relay to the destination 
is the bottleneck link. 

Remark 19: The set of power policies, B^^^ and B^^'^\ are defined by the conditions in (|93])- (|98] ). Note 
that these conditions are in general not the same as those for DF. Thus, the set B'^^^^ for the cut-set outer 
bound will in general not be exactly the same as that for the inner DF bound. However, when case 36 
maximizes the cut-set outer bounds the optimal DF = p(°^)(F) = P^^^'^{H) belongs to B^^^^ for 

both bounds. 

B. Outer Bounds Rate Region: Optimal Policy and Capacity Region 

One can similarly write the rate expressions and the KKT conditions for every point on the boundary 
of TZoB- Such an analysis will be similar to that for the K-user orthogonal MARC under DF developed 
in Section [VI-BI From Theorem [6l every point YlkeK l^^Rk on T^OB results from an intersection of two 
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polymatroids. For those cases in which the intersection is an inactive case, both the SIMO cut-set bound 
at the relay and destination and the cooperative cut-set bound at the destination are involved, and thus, 
one cannot achieve capacity. This is also true for the boundary cases. For cases 3a, 36, and 3c, in which 
the polymatroid intersection also has 2^ — 1 constraints, and hence, Kl corner points on the dominant 
K-user sum-rate face, J2keK l^kRk is maximized by a comer point of the resulting polytope. Since any 
polytope that results from some or all of the SIMO bounds will be larger than the corresponding DF inner 
bounds, the cut-set bounds are tight only when 7^2 /"2)) C TZi (E^°^\H,fj,i, fi2)^ where 

,/ii,/i2) denotes the power policy maximizing J2keicf^kRk- We summarize this observation in 
the following theorem. 

Theorem 8: The capacity region Cq-marc of an ergodic orthogonal Gaussian MARC is achieved by 
DF when for every point J2ke!cl^kRk on TZqb achieved by P^°''^(H,ij,^ , w?), 

7^2 (p(°'')(^,^i,/i2)) c 7^l (p(°')(^,/xi,/i2)) (100) 

such that P^°^\H,ixi,ii2)=P}-^^\H,^ii,^i2). Thus, Cq-marc is given by 

Co-MARC = 7^2 (Z^"'^(^,W,/i2)) = Tld [P^^'\H,fli,fl2)) . (101) 

Proposition 6 ( [3, Theorem 9]): For the case in which H_ has uniform phase fading and the channel 
state information is not known at the transmitters such that Pj^^\i£} = P^^\h) = Pk, for all k ^ T, 
Theorem H] yields the capacity region of an ergodic phase fading orthogonal Gaussian MARC as developed 
in [3, Theorem 9]. 

C. Illustration of Results 

We present numerical results for a two-user orthogonal MARC with Rayleigh fading links. We model 
the channel fading gains between receiver m and transmitter k, for all k and m, as 

H^^k = (102) 

where dm,k is the distance between the transmitter and receiver, 7 is the path-loss exponent, and Am,k 
is a circularly symmetric complex Gaussian random variable with zero mean and unit variance such that 
\Hm,k\'^ is Rayleigh distributed with zero mean and variance ^/<X^f^- For the purpose of our illustration, 
we set 7 = 3. 

Towards illustrating the sum-capacity result, we consider a two-user geometry shown in Fig. [6] For 
this geometry, in Fig. |7] we plot the inner (DF) and outer cutset bounds on the sum-rate for = 1/2 as a 
function of the relay position along the x-axis. As a result of the symmetric geometry, for every choice 
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Fig. 6. A symmetric two-user MARC geometry. 

of the relay position, both the inner and outer bounds on the sum-rate are maximized by one of cases 
3a, 36, or 3c. For each case, we use an iterative algorithm, as described in the Appendix, to compute 
the sum-rate maximizing user policies. For cases 3a and 3b, the iterative algorithm simplifies to the 
iterative waterfilling algorithm developed in [37] in which at each step the algorithm finds the single-user 
waterfilling policy for each user while regarding the signals from the other user as noise. For case 3c, 
the optimal policy at each step is still obtained by regarding the signals from the other user as noise; 
however, the user policy at each step is no longer a waterfilling solution. Finally, the optimality of DF 
when the sources are clustered relatively closer to the relay than to the destination is amply demonstrated 
in Fig . |7] The inner and outer bounds are also compared with the sum-capacity of the fading multiaccess 
channel without a relay and 6 = 1, shown by the dashed line that is a constant independent of the relay 
position. Also shown in Fig |7] are the ranges of relay positions for cases 3a, 36, and 3c for both DF and 
the cutset bounds. 

VIII. Concluding Remarks 

We have developed the maximum DF sum-rate and the sum-rate optimal power policies for an ergodic 
fading K-user half-duplex Gaussian MARC. The MARC is an example of a multi-terminal network 
for which the multi-dimensionality of the policy set, the signal space, and the network topology space 
contribute to the complexity of developing capacity results resulting in few, if any, design rules for 
real-world communication networks. For a DF relay, the polymatroid intersection lemma allowed us to 
simplify the otherwise complicated analysis of developing the DF sum-rate optimal power policies for 
the two-user and K-user orthogonal MARC and the i^-user outer bounds. The lemma allowed us to 
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Fig. 7. Plot of inner and outer bounds on the sum-rate vs. the relay position 



develop a broad topological classification of fading MARCs into one of following three types: 

i) partially clustered MARCs where a subset of all users form a cluster with the relay while the 
complementary subset of users form a cluster with the destination 

ii) clustered MARCs comprised of either sources-relay or relay-destination clustered networks, and 

iii) arbitrarily-clustered MARCs that are a combination of either the two clustered models or of a 
clustered and a partially clustered model. 

The optimal policies for the inner DF and the outer cutset bounds for the orthogonal half-duplex 
MARC model studied here lead to the following observations: 

• that DF achieves the sum-capacity of a class of source-relay clustered orthogonal MARCs for which 
the combined link from all sources and the relay to the destination, i.e., the link achieving the K-usei 
sum-rate at the destination, is the bottle-neck link. Furthermore, DF achieves the capacity region 
when for every weighted sum of user rates, the limiting bound is the weighted rate-sum achieved 
at the destination. 

• that for this sum-capacity achieving case, the optimal user policies for both the orthogonal and non- 
orthogonal half-duplex MARCs are multi-user opportunistic waterfilling solutions over their links 
to the destination and the optimal relay policy is a water-filling solution over its direct link to the 
destination. 
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• and that for the remaining classes of MARCs, the optimal users policies are waterfilling and non- 
waterfiUing solutions for the partially clustered and arbitrarily clustered models, respectively. 

For the partially clustered cases, we have showed that the optimal pohcy for each user is multiuser 
waterfilling over its bottle-neck Unk to one of the receivers. Thus, the users that are clustered with the 
destination are forced to transmit at a lower rate to allow decoding of their signals at the relatively distant 
relay. Our results suggest that a useful practical strategy for the partially clustered topologies may be 
to allow those distant users that present little interference at the relay to communicate directly with the 
destination. 

The optimal relaying strategy for all except the capacity achieving clustered case described above 
remains open. Given the complexity of finding the optimal signaling schemes for a given performance 
metric in multi- terminal networks, a natural extension to this work could be to understand the gap in 
spectral efficiency between DF and the cutset outer bounds for fading MARCs. Such bounds have been 
developed recently for time-invariant interference chaimels and relay chaimels in [38] and [39], respec- 
tively, and for fading Gaussian broadcast chaimels with no chaimel state information at the transmitter 
in [40]. 

Our analysis can also be extended to study the more general orthogonal half-duplex MARC model 
where the sources transmit on both orthogonal channels while the half-duplex relay is limited to receiving 
on one and transmitting on the other. The half-duplex relay receives signals from the sources on one 
of the bands while the destination receives it in both bands. For the special case where the destination 
can only receive in the band used by the relay to transmit, we obtain a multiple-access version of the 
orthogonal relay channel studied in [12]. Irrespective of the receiving capabihties of the destination, for 
this more general orthogonal model, each source transmits two signals, one for each band, subject to an 
average power constraint over both bands. 

Thus, for this general model, as in [12], one can consider a more general decoding scheme of 
partial decode-and-forward (PDF) where each source transmits two independent messages, one on each 
orthogonal channel (see also [10]). As a result, the analysis does not simplify to studying an intersection 
of two polymatroids as it does for DF. However, analogous to the time-invariant (non-fading) case, we 
expect that PDF will simplify to the sum-capacity optimal DF for the special case in which all the 
sources and the relay are clustered. While this general orthogonal model is useful to study for the sake 
of completeness, the model we study here abstracts practical multi-hopping architectures and provides 
insights into network architectures and topologies where using a decode-and-forward relay is beneficial. 

Finally, a note on complexity: our theoretic analysis distinguishes between all possible polymatroid 
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intersection cases in determining the optimal policy for a K-user system and therefore has a complexity 
that grows exponentially in the number of users. In practice, however, for two intersecting polymatroids 
the maximum of a weighted sum of rates and the optimizing policies can be computed using strongly 
polynomial-time algorithms [26, Theorem 47.4]. 



Appendix 
Proof of Theorem [3] 

The sum-rate maximizing DF power policy P* {H_) in Theorem [3] is obtained by sequentially deter- 
mining the power policies P^*) (H) and pC'") (^) that maximize the sum-rate for cases i and {l,n), 
respectively, over all P [H_) € V, until one of them satisfies the conditions for its case. We consider each 
case separately starting with case 1. 

Case 1: This case occurs when the power policy P_{H_) € Bi is such that the intersection of the relay 
and destination rate regions belongs to the set of inactive cases (see Fig |2l). The Lagrangian for this 
sum-rate maximization is given by 

= SW - Y^^^^ i.,E [Pfc {H) - Pfc] + Y,^^^ ^kPk (H) (103) 

where, for all k ^ T, Uk are the dual variables associated with the power constraints in ([5]), > are 
the dual variables associated with the positivity constraints P^. {H_) > 0, and 



'ec (\Hd,i\^ Pi {H)/e) + ec (\Hd,r\^ p, m/e)] + E\ec (\Hr,2\^ m/e 



(104) 



The optimal policy P*^^-* {W) maximizes (11031) if it belongs to the open set Bi defined by the conditions 

R{i}4 {P^'^ (H)) < Rftlr {P^'^ (H)) and P{2},r {P^'^ (K)) < Pgf , (p(^) {H)) (105) 



where 



^v"" {P{K)) = OHXi; Yr\H) = E 



ec 



RfX^nH)) = ei{X2;Ya\H)- 
The KKT conditions for (11031 ) simplify to 



E 



ec 



\Hr,l\^ Pi{H)/e 



dC 
dPk{h) 



fl^^ - i^fc In 2 < 0, with equality for P^ (h) >0, k = l,2,r 



(106) 
(107) 

(108) 
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where 



/^'^ = (TT^frPl^ (fc,m) = (l,d),(2,r), (109) 

It is straightforward to verify that these KKT conditions result in 

(/^) = (i;^-ct)^ (A:,m) = (l,d),(2,r) (111) 



and 



Ca^e 2: With 1^,^. and as the dual variables associated with the power and positivity constraints on 
Pk, respectively, the Lagrangian for this case is 

£(2) = 5(2) _ ^^^^ ^^JE [p^ (^) _ p^] + ^^^^ A.Pfc (^) (113) 

where 

S^'^ = R{l},r (£ (K)) + R{2}4 (£ (£)) (1 14) 



E 



(115) 



ec Pi {H)/9^ + 0c {\Hd,2\'' P2 m/e) + (\Hd,r? Pr m/'i 

The optimal policy 

p(2) (p;) 

maximizes (11031 ) if it belongs to the open set 132 given by the conditions 

R{i}APiH)) < Rflf^APm and i?{2M(£(£))<i?gf,,, (£(£)) (116) 

where R'^2}r ^^'^ ^{I'f d given by (11061 ) and (11071 ). respectively, after replacing the user indices 1 by 
2 and 2 by 1. Note that 5(2) and £(2) 

can be obtained from S^-^^ and £.^^\ respectively, by interchanging 
the user indices. Thus, the optimal pj:^^ (H) and pP (H) are given by (lllll ) and (11121 ). respectively, 
with {k,m) = (l,r), (2, d) provided P^^) (H) satisfies (fTT6l) . 

Cai'e i.' Consider the three cases 3a, 36, and 3c shown in Fig. [3l The sum-rate optimization for all 
three cases is given by 

ma.xmm{Ric,r,Rjc,d) (117) 

subject to average power and positivity constraints on P^, for all k. Recall that we write Rfcj to denote 
the sum-rate bound at receiver j where the two bounds at the relay and destination are given by (ITSl) and 
([T9l ). respectively. We write to denote the open set consisting of all P {H_) G V that do not satisfy 
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(11051) and (11161 ) either as strict inequalities, i.e., do not satisfy the conditions for cases 1 and 2, or as 
a mixture of equalities and inequalities, where by a mixture we mean that a subset of the inequalities 
in ( 11051 ) and (11 161 ) are satisfied with equality. We will later show that such sets of mixed equalities 
and inequalities in (11051 ) and (II 161 ) corresponds to conditions for the various boundary cases (see also 
Figs. |4] and [5]). Thus, P {H_) G S3 only when it does not satisfy the conditions for the inactive and the 
active-inactive boundary cases. By definition, B3 = B^a U B^^, U ^sc, where Bi, i = 3a, 3b, 3c, is defined 
for case i below. 

The optimization in (II 17b is a multiuser generalization of the single-user max-min problem studied in [6] 
(see also [22, Sec. 3.1]) for the orthogonal single-user relay channel. In [6], the authors use a technique 
similar to the minimax detection rule in the two hypothesis testing problem (see for e.g., [41, II.C]) to 
show that the max-min problem simplifies to optimizing three disjoint cases in which the maximum rate 
is achieved either at the relay or at the destination or at both (boundary case). The classical results on 
minimax optimization also applies to the multi-user sum-rate optimization in (11171 ). and thus, the optimal 
policy P(') (H), i = 3a, 3b, 3c, satisfies one of following three conditions 

Case 3a: Rjc,r\p(3-)(^H_) < Rk4\p^^'^){h_) 



Case 3b: R}c,r\pt.^'')(H_) > RK,d\pi^''>(H_) 



(118) 
(119) 

Case 3c: i?/c,r|p(3<=)(_H;) = RK,d\pi^'=){H_)- (120) 

Note that the conditions in (II 18I )- (I120I ). evaluated at any P G B3, are also conditions defining the sets 
fSsa, S^b, and B-^c, respectively. Before detailing the optimal solution for each of the above three cases, 
we write the Lagrangian for case i as 

~^ 3a, 36, 3c (121) 

XkPk m > (122) 



k£T k&T 



where Vk and > are dual variables associated with the average power and positivity constraints on 
Pk, respectively, 

E 



RK,r 



E 



C E \Hr,k\' Pkim/i 

\k=l 

cf E \Hd,k\^ Pk{H)/i 

\k=l 



(1 - a) Ric,r + aRicd 



3a 

36 
3c, 



(123) 



and a is the dual variable associated with the equality (boundary) condition in (11201 ). The resulting KKT 
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conditions are given by 



9^ - F^'^ = f^^ - Ukln2 < 0, k = l,2,r,i = 3a,3b,3c (124) 



dPk{h) k Jk 

where (11241 ) holds with equaUty for Pk {H_) > and for A; = 1, 2 



Jk 



Kkf / {i+ E \K,k\^ Pk{h)/e] i = 3a 

k=l 



\hd,kf ^{^1 + pjhd,kf Pkih)/9^ i = 3b (125) 
I (l-a)/f)+a/f) i = 3c. 



and /i'^"-' = /J^''^ = fr^\ fr^^^ = afr^\ We now present the optimal policies and sum-rates for each 
case in detail. 

Case 3a: For this case, the KKT conditions in (|124l) and (11251 ) depend only the sum-rate and channels 
gains of the two users at the relay. Thus, the problem simplifies to that for a MAC channel at the relay 
and the classic multiuser waterfilling solution developed in [27], [28] applies. From (1124b . the optimal 
user policies are 

With the exception of the equality condition in (11261 ). the optimal policies are unique, i.e., the optimal 
p^3a) ^j^^ ^^^^ ^ (11261 ) is an opportunistic water-filling solution that exploits the fading diversity 
in a multiaccess channel from the sources to the relay. If the channel gains are jointly distributed with a 
continuous density, the equality condition occurs with probability 0. Furthermore, even if the distributions 
were not continuous, one can choose to schedule one user or the other when the equality condition is 
met, thereby maintaining the opportunistic allocation poUcy. Finally, the optimal power policy at the relay 
is not explicitly obtained from C^^^ in (11211 ) as for this case S^^^ is the sum-rate achieved by the sources 
at the relay. However, since the sum-rate at the relay for this case is smaller than that at the destination, 
choosing the optimal waterfilling policy at the relay that maximizes the relay-destination link preserves 
the condition for this case, and thus, pj:^"^ (H) IS given by (fml) . When P^^a) ^jf^ ^ requirement 
of satisfying (1118b . i.e., (H_) g Bsa, simpUfies to a threshold condition Pr > Pu (^1,^2) where 

Pk, A; G T, is defined in ([5]) and the threshold Pu [Pi,P2} is obtained by setting (II 181 ) to an equality. 
When P(3a) g but P^^") (H) ^ Bsa, Ri + R2 is maximized by either case 3b or case 3c. For 
as argued in Section HVl the sum-rate is not maximized by any P_{H_) G B3. 
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Case 3b : The optimal policy Pjf''^ {H_) at user k for this case satisfies the KKT conditions in (11241 ) with 
ff = ff^^ in dmll. As with case 3a, here too, the optimal policy is an opportunistic water-filling solution 
and is given by (11261 ) with the subscript V changed to 'd' for all k and with the superscript i = 3b. 
Further, for the relay node, the optimal P^^^^ (H) satisfies the KKT conditions in (11081 ). i.e., fi^^'^ = fi^\ 
and is given by the water-filling solution in (II 121 ). Finally, for p(^^) (^H_) G B3, the requirement P^^**) (H) 
G B^h simplifies to satisfying the threshold condition < Pi (-Pi, -P2) where Pi (Pi, P2) is determined 
by setting ( 11191 ) to an equality. 

Case 3c (equal-rate policy): The optimal policy Pjf^^ {H_) at user k for this case satisfies the KKT 
conditions in (11241) for /^*^ = fjf'^^ in ( 11251 ). The function fj^'^^ in ( 11251 ) is a weighted sum of fjf'^^ and 
fj^^^ where the Lagrange multiplier a accounts for the boundary condition in (11201 ). Substituting fj^^^^ 
in ( |125b in (11241 ). we have the following KKT conditions 



a\hr.M' _^ {i-a)\h^.,\' <Ukln2 with equality for Pk (h) > 0, k = 1,2 



(127) 



k = l 



which impUes 



/fV^i>/fV^2 Pr^(/i) = frootofpf'=)|p,=o) ,PtHh) 



/f 'V^i < ft^H Pf (h) = 0, Pt^ (h) = (root of |p,=o) (128) 
/f'V^i = ft^H Pt^ (h) and Pt^ (h) satisfy /f = ^,ln2 
where Pjf^^ is defined in (11241 ). Determining the optimal Pjf'^^ (h), k = 1,2, requires verifying each one 
of the three conditions in (11281 ). Note that in contrast to case 3a (and case 35 with 'r' replaced in (11261 ) 
by 'd'), the opportunistic scheduling in (11281 ) also depends on the user policies in addition to the channel 
states. Furthermore, the optimal solutions Pjf^^ {H_) do not take a water-filling form. Thus, for a given 
Pi (h), P2 (h) is given by 

P2 {h) = positive root x of (11301 ) if it exists, otherwise (129) 

where the root x is determined by the following equation: 

+ . =.2ln2. (130) 



l + \llr,l\ — 5-^ + |'ir,2| l + \lld,k\ —g-^ + \n'd,k\ 

Using P2 (h) given by (11301) . Pi (h) is obtained as the root of 

a\hr,i\'^ (1 - a) IViP = ,,]n9 n3n 

Thus, for all h, starting with an initial Pi (h), we iteratively obtain Pi (h) and P2 (h) until they converge 
to Pf (H) and P^^"^ (H). The proof of convergence is detailed below. Finally, the optimal policies are 
determined over all a G [0, 1] to find an a* that satisfies the equal rate condition in (1120b . 
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Proof of Convergence: The proof follows along the same lines as that detailed in [22, p. 3440] and 
relies on the fact that the maximizing function S^^'''^ in (fT23]) is a strictly concave function of Pi (H) 
and P2 {H_) and is bounded from above because of the power constraints at the source and relay nodes. 
In each iteration, the optimal Pi {H_) and P2 {H_) are the KKT solutions that maximize the objective 
function. Thus, after each iteration, the objective function either increases or remains the same. It is easy 
to check that for a given Pi {H_) the objective function is a strictly concave function of P2 {H_), and 
thus, (11301 ) yields a unique value of P2 {H_). Furthermore, the objective function is also a strictly concave 
function of Pi {H_) for a fixed P2 {H_). Thus, as the objective function converges, (Pi {H_) , P2 {H_)) also 
converge. Finally, Pi {H_) and P2 {H) converge to the solutions of the KKT conditions, which is sufficient 
for (Pi (H) , P2 (H)) to be optimal since the objective function is concave over all P_{H_) G V. 

Finally, since fj^^^ = afr^\ the relay's optimal policy simplifies to the water-filling solution given by 

^"'*^'^(^-^)'- 

Case 4: {Boundary Cases): Recall that we define the sets Pj, i = 1,2, 3a, 36, 3c, as open sets to ensure 
that an optimal P* maximizes the sum-rate for a case only if it satisfies the conditions for that case. 
Since an optimal policy can lie on the boundary of any two such cases, we also consider six additional 
cases that lie at the boundary of an inactive and an active case. These boundary cases result when the 
conditions for an inactive case 1,1 = 1, 2, and an active case n, n = 3a, 36, 3c, are such that the sum-rate 
is the same for both cases. We consider each of the six boundary cases separately and develop the optimal 
p{i,n) for each case. The requirement that the optimal pC'") (H) satisfies the condition S"^'^ = 5^"^ 



for the boundary case (/,n) simplifies to 

case (1, 3a) P{i},d + R{2},r = RK.,r < RK.,d (133) 

case (1, 36) R^^^d + R{2},r = Ric,d < Ric,r (134) 

case (1, 3c) P{i}^rf + R{2},r = Rjc,r = Rka (135) 

case (2, 3a) R{i}^r + R{2},d = RK.,r < Rk.4 (136) 

case (2, 36) P|i},,, + R{2},d = Rjc,d < Rjc,r (137) 

case (2, 3c) P{i}^^ + R{2},d = Ric,d = RjCr (138) 



where the conditions in (I133b - (ll38b are evaluated at the appropriate p(''") [H_). Note that the conditions in 
(|133I )- (I138I ) also define the conditions for the sets ^(i,3a) through i3(2,3c), respectively. Using (|133I )- (I138I ). 
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we write the Lagrangian for all boundary cases except cases (1,3c) and (2,3c) as 

^(z,n) ^ ^gii) + (!_„) 5{n) _ ^ [Pk {H) -Pk]+Y: ^kPk {H) / = 1, 2, n = 3a, 36 

fcGT fcer 

(139) 

AfcPfc {H) > (140) 
and the Lagrangian for cases (1,3c) and (2,3c) as 

jriW ^ ^ ^^g{3a) + (1 _ „^ _ „2) 5(3^) - [Pk m - Pk] 

ker 

+ E ^kPk{H), 1 = 1,2 (141) 
XkPk m > (142) 

where Uk and Afc > are dual variables associated with the average power and positivity constraints on 
Pk, respectively. The variable a is the dual variable associated with all boundary cases with a single 
boundary condition while ai and a2 are the dual variables associated with cases (1,3c) and (2,3c). The 
resulting KKT conditions, one for each Pk (h), k = 1,2, r, are 

Case (/, n ^ 3c) : |g£ = fj!'^^ = aff + {1 - a) ff^ < z.. In 2 (143) 

Case (/, n = 3c) : fgg = /(''"^ = ai/f + aa/f + {1 - ai - a^) /f <Vk\n2 (144) 

where /^'^ and /^"^ are as defined earlier for cases I and n and equality holds in (11431 ) and (11441 ) for 
Pk {bd > 0, for all h. We now present the optimal policies for each case separately. 
Case (1,3a): From (11431 ). the KKT conditions for this case are 



p(l,3a) _ a|fed.r| 
r 

which implies 



= ^ 2 with equaUty if P, (h) > (147) 



iTl^f^ A(/i) = frootof Ff'''^V.=o) ,^2(Zi) = 



1/1 ' 1^2 

r(l,3a) r(l,3a) 



^ < ^ A (A) = 0, P2 (h) = (root of F^'''^^ |p,=o ) (148) 



1/2 

P(l,3a) r(l,3a) r(1.3a) r(l,3a) 

^ = ^ ^1 (^) and ^2 (/i) satisfy A_ = L_ 



where 



^{i,n) ^ ^(i,n) _ In 2 < 0, for all {l,n) . (149) 
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As in case 3c, the optimal policies take an opportunistic non-waterfilling form and in fact can be obtained 
by the iterative algorithm described for that case. Finally, from (fT47]) . the optimal pj^'^"^ (H) IS given 
by (im . 

Case (1,36): The analysis for this case mirrors that for case (1,3a) and the optimal user policies are 
opportunistic non-water-filling solution given by (11481 ) with fj^^^ replaced by fj?^\ k = 1,2. On the 
other hand in contrast to case (1, 3a) where fi^'^^ = 0, since both fP and /i^''^ are non-zero, the optimal 
relay policy P^^^'^"^ = 

Case (1,3c): For this case, the KKT conditions in (11441 ) involves a weighted sum of fj:\ fj?"'\ and 
/f^^ Thus, for k = l,2, {k,m) = (l,(i),(2, r), we have the KKT conditions 

= aiff^ + as/f + (1 - ai - aa) /f < ui In 2 with equality if Pi {h) > (150) 

f^^'^^^ = + 02 /f"^ + (1 - /?) (1 - a) /2^^^^ < In 2 with equality if P2 (/i) > (151) 

= (1 - 02) IVrP/c (^IV^pPr/^) <z^rln2 with equality if P^ (/i) > (152) 

where ai and 02 are the dual variables associated with the equalities Ric,d = R{i},d + ^{2},r ^^d 
RK,d = Ric,r, respectively, in (11351 ). From (1150b and (11511 ). one can verify that the optimal user policies 
are opportunistic non-waterfilling solutions given by (11481 ) with the superscript (1, 3a) replaced by (1, 3c). 
Finally, pj^'^^^ {H) is given by the water-Ming solution in (1112b with a replaced by (1 — 02). 
Case (2,3a): The optimal user policies for this case and the KKT conditions they satisfy are given by 
(11451 ). (11461 ). and (1148b when fj^^^ is replaced by /(^\ for all k, and g^^ is superscripted by (2, 3a). Thus, 
here too, the optimal user policies are opportunistic non-water-filling solutions. The optimal relay policy 

same as that obtained in case (1, 3a). 
Case (2,36): The optimal user policies p^^'^''^ (p^), A; = 1,2, are again opportunistic non-water-fiUing 
solutions and are given by (11451) . (fT46l) . and (fT48l) when ff^ and /f ''^ are replaced by and f<^''K 
respectively, for all k, and (7^^ is superscripted by (2,36). The optimal relay policy pj^''^^^ (P) is the 
same as that for case (1,36). 

Case (2,3c): The optimal policy vector p(2'3c) (^ff^ jg {j^g same as that for case (1,3c) with /^^^ is 
replaced by f^'^\ for all k, and with the superscript (2,3c). 
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